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ABSTRACT 



^ 
Ph 



^^ We investigate the origin and physical properties of Ovi absorbers at low redshift 

rj (z - 0.25) using a subset of cosmological, hydrodynamical simulations from the Overwhelm- 

ingly Large Simulations (OWLS) project. Intervening Ovi absorbers are believed to trace 
shock-heated gas in the Warm-Hot Intergalactic Medium (WHIM) and may thus play a key 
role in the search for the missing baryons in the present-day Universe. When compared to 
Q observations, the predicted distributions of the different O vi line parameters (column density 

H Novi, Doppler parameter ^ovi, rest equivalent width W^) from our simulations exhibit a lack of 

c/5 strong O vi absorbers, a discrepancy that has also been found by Oppenheimer & Dave (2009). 

^ ^^ This suggests that physical processes on sub-grid scales (e.g. turbulence) may strongly influ- 

ence the observed properties of Ovi systems. We find that the intervening Ovi absorption 
^vq arises mainly in highly metal-enriched (10"' <K ZIZq < 1) gas at typical overdensities of 

J> 1 <s: p/ {p} < 10^. One third of the Ovi absorbers in our simulation are found to trace gas 

(^ at temperatures T < 10^ K, while the rest arises in gas at higher temperatures, most of them 

■^ around T = j^q^^'^*'^ K. These temperatures are much higher than inferred by Oppenheimer & 

OO Dave (2009), probably because that work did not take the suppression of metal-line cooling 

Cn by the photo-ionising background radiation into account. While the O vi resides in a similar 

[~ — region of (p, r)-space as much of the shock-heated baryonic matter, the vast majority of this 

f^ gas has a lower metal content and does not give rise to detectable O vi absorption. As a con- 

^^ sequence of the patchy metal distribution, O vi absorbers in our simulations trace only a very 

^~~^ small fraction of the cosmic baryons (< 2 percent) and the cosmic metals. Instead, these sys- 

►^ tems presumably trace previously shock-heated, metal-rich material from galactic winds that 

• ^H is now mixing with the ambient gas and cooling. The common approach of comparing O vi 

/\^ and H i column densities to estimate the physical conditions in intervening absorbers from 

H QSO observations may be misleading, as most of the H i (and most of the gas mass) is not 

physically connected with the high-metallicity patches that give rise to the O vi absorption. 

Key words: cosmology: theory — methods: numerical — intergalactic medium — quasars: 
absorption lines — galaxies: formation 



1 INTRODUCTION the baryons at redshift z = 3 (Rauch et al. 1997; Weinberg et al. 

... . 1997; Schaye 2001). However, at z = the fraction of baryons 

Diiiuse ionised sas in the intergalactic medium (IGM) represents .,..,. ^ . , , , ^^ a^ 

. . .^ residing in the Lya forest is strongly reduced to ~ 30 - 40 per 

the maior baryon reservoir in the Universe at any redshift. From . ^r. . . i or,n..x i -i » »i .■ ^, , , i i n 

^ cent (Penton et al. 2004), while at the same time the (observable) 

observations of the Lya forest line density at high redshift in quasar 

(QSO) and active galactic nuclei (AGN) spectra it can be deduced 



amount of baryons in condensed structures (i.e., baryons in galaxies 

and galaxy clusters) has increased to only a few per cent (Fukugita 
that the photo-iomsed IGM contains more than ninety per cent of r^nnt-^ t^, i ^- .i • j- . ^i . ■ -n » r .■ 

^ ^ ^ 2004). These observations thus indicate that a signmcant fraction 

of the baryons formerly residing in the photo-ionised IGM have 

"disappeared". Cosmological simulations have predicted that, as a 

* E-mail: tepper@astro.physik.uni-potsdam.de 
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result of the large-scale structure formation in the Universe, most of 
these "missing baryons" have moved into a hot, shock-heated inter- 
galactic gas phase (Cen & Ostriker 1999; Dave et al. 2001; Bertone 
et al. 2008). This shock-heated intergalactic gas phase is referred to 
as the Warm-Hot Intergalactic Medium (WHIM) and is expected to 
have characteristic temperatures in the range T ~ 10' - 10 '' K and 
densities of «h ~ lO"'' - lO^^cm"^. Constraining the distribution 
and physical properties of the WHIM is important to understand 
how the gas (and the metals) are transported from galaxies into the 
IGM and recycled into galaxies, and what role the shock-heated 
IGM has for the evolution of galaxies in the local Universe. 

Since emission from such a thin plasma is extremely dim (Fur- 
lanetto et al. 2004; Bertone et al. 2009, 2010) and the hydrogen 
is almost fully ionised, the analysis of highly ionised heavy ele- 
ments (in particular the high ions of oxygen, O vi O vii , and O viii) 
in Ultra- Violet (UV) and X-ray absorption against distant extra- 
galactic background sources has become the leading method to 
study the properties and baryon content of the WHIM at low red- 
shift (for a recent review see Richter et al. 2008). Most promising is 
the search for five-times ionised oxygen (O vi) in the UV spectra of 
low-redshift AGN, as obtained with space-based UV spectrographs 
such as Hubble Space Telescope (HST) Space Telescope Imaging 
Spectrograph (STIS) and the Far Ultraviolet Spectroscopic Ex- 
plorer (FUSE) {e.g. Tripp et al. 2000; Richter et al. 2004). Oxygen 
is a relatively abundant element with two strong O vi transitions at 
/il031.9A and /il037.6A, and so far, more than 50 intervening O vi 
absorbers at low redshift have been identified and analysed (e.g. 
Tripp et al. 2008; Danforth & ShuU 2008). 

In spite of the large O vi samples obtained to date, there is still 
no general consensus about the physical conditions of the gas giv- 
ing rise to O vi absorption. While Danforth & Shull (2008) find that 
O VI (and associated N v) are reliable tracers of coUisionally ionised 
gas at temperatures 10' K < T < 10* K (i.e., the low-temperature 
WHIM), Thom & Chen (2008a) argue that Ovi arises mainly in 
photo-ionised gas at temperatures T < 10' K. Similarly, Tripp et al. 
(2008) find that well-aligned Ovi-Hi absorbers have line widths 
that are consistent with photo-ionised gas. Nevertheless, these au- 
thors show that more than half of their O vi absorbers are complex, 
i.e., multi-phase, and could thus trace both cold photo-ionised gas 
and lower metallicity coUisionally ionised gas at T > lO' K. 

The interpretation of the abundance and nature of interven- 
ing Ovi absorbers arising in coUisionally ionised gas in terms of 
the distribution and baryon content of the WHIM is not straight- 
forward. In coUisional ionisation equilibrium (CIE; as usually as- 
sumed for a shock-heated plasma like the WHIM), O vi predomin- 
antly traces gas at T ~ 3 x lO' K, while in the higher temperature 
regime of the WHIM oxygen is further ionised to O vn and O viii, 
observable only in the X-ray band for which high-quality spectral 
data are sparse. Thus, only a (small) fraction of the WHIM can 
in principle be traced by intervening O vi absorbers. In addition, 
several authors {e.g. Gnat & Sternberg 2007) have argued that the 
assumption of CIE is not valid in case of the WHIM {i.e., the gas 
may be out of an ionisation equilibrium) and more recently, Wi- 
ersma et al. (2009a) have shown that photo-ionisation of the WHIM 
strongly influences the cooling efficiency of the gas in the temper- 
ature range of interest. These physical arguments, together with the 
not well established oxygen abundance of the WHIM at low red- 
shift, indicate that a reliable estimate of the baryon budget of the 
WHIM (i.e., the amount of ionised hydrogen) from Ovi observa- 
tions requires a deep understanding of the physical conditions of 
the gas. 

In this paper, we make use of a set of numerical simula- 



tions from the Overwhelmingly Large Simulations (OWLS) pro- 
ject (Schaye et al. 2010) to study the physical properties and baryon 
content of intervening O vi absorbers at low redshift. The main ad- 
vantage of these simulations compared to previous ones is the im- 
plementation of important physical processes that have been largely 
ignored in earlier studies. As will be shown here, the influence of 
the photo-ionisation on the cooling function of the WHIM (Wi- 
ersma et al. 2009a) represents a particularly important aspect for 
a correct interpretation of the observed properties of the O vi ab- 
sorbers and their role for our understanding of the WHIM. 

Even though we do not attempt to tune our simulations to re- 
produce any Ovi observables in this paper, we will compare our 
results to observations where possible. Moreover, a thorough ana- 
lysis of variations with respect to our fiducial model and the effect 
of post-run physics variations on the resulting physical properties 
of the simulated O vi absorbers are left for a future paper. Note that 
we will often refer to Oppenheimer & Dave (2009) and make com- 
parisons between our and their results, since Oppenheimer & Dave 
(2009) is the most comprehensive study on O vi absorbers in simu- 
lations to date (see also Cen & Chisari 2010). As we will show, the 
fact that Oppenheimer & Dave (2009) assumed CIE when comput- 
ing the contribution of heavy elements to the radiative cooling rates 
will likely have affected their conclusions. 

This paper is organised as follows: in Sec. 2 we briefly intro- 
duce the simulation we use, i.e., our fiducial model. We describe 
how we create our synthetic spectra and present their analysis in 
Sec. 3. In Sec. 4 we discuss a comparison of our results to ob- 
servations, while in Sec. 5 we analyse the physical properties of 
the O VI gas in our simulation in detail. We present and discuss our 
conclusions in Sec. 6. A detailed analysis of the effect of the signal- 
to-noise ratio, simulation box size, and resolution on the O vi line 
parameter distributions are left for the Appendix. 



2 SIMULATIONS 

The simulations used in this work are part of a large set of 
cosmological simulations that together comprise the OWLS 
project, described in detail in Schaye et al. (2010, and references 
therein). Briefly, the simulations were performed with a signi- 
ficantly extended version of the A'-Body, Tree-PM, Smoothed 
Particle Hydrodynamics (SPH) code gadget iii - which is a 
modified version of gadget n (last described in Springel 2005) 
- a Lagrangian code used to calculate gravitational and hydro- 
dynamical forces on a system of particles. The initial conditions 
were generated from an initial glass-like state (White 1996) with 
CMBFAST (version 4.1; Seljak & Zaldarriaga 1996) and evolved 
to redshift z = 127 using the Zel'Dovich (1970) approximation. 
The reference model of OWLS adopts a flat ACDM cosmology 
characterised by the set of parameters {fim, fib^ f^^, cg, n^, h} = 
{0.238,0.0418,0.762,0.74,0.95,0.73) as derived from the 
Wilkinson Microwave Anisotropy Probe (WMAP) 3-year data' 
(Spergel et al. 2007). The reference model includes prescriptions 
for star formation (SF) and for feedback from core collapse 
supernovae (SNe) described in Schaye & Dalla Vecchia (2008) and 
Dalla Vecchia & Schaye (2008), respectively. The implementations 
of radiative cooling and heating are described in Wiersma et al. 

' These parameter values are largely consistent with the WMAP 7-year 
results (Jarosik et al. 2010), the largest difference being the value of erg, 
which is 2 o" lower in the WMAP 3-year data than allowed by the WMAP 
7-year data. 
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Table 1. List of simulations -which assume the reference model of the OWLS project- used in this study. Our fiducial run REF_L050N512 (shown in bold) is 
used for the comparison to observations and our general analysis, whilst all other simulations are used in the Appendix when addressing the convergence of our 
results with respect to box size and resolution. The columns show": simulation identifier''; comoving box size (L); number of dark matter particles {N; there 
are equally many baryonic particles); baryonic particle mass (nn,); dark matter particle mass (mdm); comoving (Plummer-equivalent) gravitational softening 
(fcom); maximum physical softening (fprop); final redshift (zend)- 



Simulation 


L 


N 


nib 


'«dm 




^com 




^prop 


■Send 




(h-^ Mpc) 




(A-i Mo) 


(A-1 Mo) 


(h- 


'kpc) 


(h- 


'kpc) 




REFX025N128 


25.00 


128^ 


8.7 X 10' 


4.1 X 10* 




7.81 




2.00 





REFX050N128 


50.00 


128^ 


6.9 X 10* 


3.3 X 10' 




15.62 




4.00 





REFJM50N256 


50.00 


2563 


8.7 X 10' 


4.1 X 10* 




7.81 




2.00 





REFX050N512 


50.00 


5123 


1.1 X 10' 


5.1 X 10' 




3.91 




1.00 





REFJ.100N512 


100.00 


5123 


8.7 X 10' 


4.1 X 10* 




7.81 




2.00 






" For a complete list of the OWLS runs as well as the detailed description of the corresponding model see Schaye et al. (2010). 

'' The name convention is [model] JjcxxNyyy, where 'xxx' and 'yyy' are, respectively, the box size in comoving A"' Mpc and cube root of the number of 

particles per species (dark matter or baryonic). 



(2009a) and summarised below (cf. Sec. 5.3). For a thorough de- 
scription of the implementation of stellar evolution in the reference 
model we kindly refer the reader to Wiersma et al. (2009b). 

We analyse the physical properties of O vi absorbers identi- 
fied in simulated spectra, focusing on a high resolution simula- 
tion run which assumes the reference model of the OWLS dubbed 
REF_L050N512. This simulation was run down to z = in a peri- 
odic box of size L = 50 comoving h^ ' Mpc, using 5 1 2 ' dark matter 
particles and equally many baryonic particles. This run has been 
shown to have a large enough box size and high enough resolution 
to obtain a converged prediction for the cosmic star formation his- 
tory for z < 4, which, however, drops off less rapidly with time than 
is observed for z < 1 (Schaye et al. 2010). Also, the predicted metal 
mass distribution for different components (stars, star-forming gas, 
non-star-forming gas) and gas phases from this particular run is 
shown to be converged for z < 2 and agrees broadly with observa- 
tions (Wiersma et al. 2009b, see discussion in their Appendix C). 
We refer the reader to Schaye et al. (2010, and references therein) 
for a more detailed description of our fiducial run (and variations 
thereof). In Table 1 we summarise its basic numerical properties, 
together with the corresponding properties of other runs assum- 
ing the reference model. These runs are discussed in Appendix B, 
where we address the convergence of our fiducial run with respect 
to box size and resolution, using various quantities derived from 
synthetic spectra as explained in the following section. 



3 SPECTRAL ANALYSIS 

In this section we briefly describe our method to compute synthetic 
spectra using the spectra generating package specwizard written by 
Schaye, Booth, & Theuns, which follows the approach described 
in Theuns et al. (1998, their Appendix A4). A detailed analysis of 
the effect of the signal-to-noise ratio (S/N) on the distributions of 
line parameters (column density, Doppler parameter, rest equival- 
ent width) derived from the spectra is presented in Appendix A, and 
the main results are only summarised here for brevity (see below). 



a given point on opposite faces of the simulation box, and the col- 
lection of SPH particles with projected distances to this line smaller 
than their corresponding smoothing length. The next step is to cal- 
culate the ionisation balance for each SPH particle as a function of 
redshift, density, and temperature, which we do using precomputed 
tables obtained with the photoionisation package cloudy (version 
07.02 of the code last described by Ferland et al. 1998), assuming 
the gas is exposed to the Haardt & Madau (2001) model for the 
X-Ray/UV background radiation from galaxies and quasars. 

A synthetic spectrum is then computed as follows. For a 
chosen pixel size Ax (in proper length units), the sightline is di- 
vided into Wpix = [a(z) L/h] /Ax bins in distance, where h and 
a(z) are the normalised Hubble constant and the expansion factor 
at the box's redshift z, respectively. Given the positions, peculiar 
velocities, and temperatures of all relevant SPH particles as well as 
the densities for each ion species (e.g., Ovi), we compute the ion 
density-weighted temperature and peculiar velocity for each spe- 
cies at each bin according to the SPH interpolation scheme. Proper 
distance bins Ax along the sightline are transformed into velocity 
bins Av via Ai; = H(z) Ax, ion number densities «!„„ into column 
densities Nio„ = ni„„ Ax, and temperatures into Doppler widths. The 
optical depth t(v) at each pixel for a particular transition is com- 
puted assuming a thermal (i.e., Gaussian) profile, taking peculiar 
velocities into account. Finally, the optical depth spectrum is trans- 
formed into a continuum-normalised flux via F(v) = exp[-T(j))]. 

We generate 1000 random sightlines through the simulation 
box of our fiducial run at 5 different redshifts spanning the range^ 
z = 0.0 - 0.5, from which we synthesise a total of 5000 continuum- 
normalised spectra containing absorption by the strongest trans- 
ition ^1031.93 A of the O vi doublet (^1031.93 A, ^1037.62 A). 
In order to mimic (and thus compare our results to) observations 
performed by HST/STIS as closely as possible, we convolve our 
spectra with a instrumental Gaussian Line-Spread Function (LSF) 
with a Fufl-Width-At-Half -Maximum FWHM = 7 km s~' , and res- 
ample our spectra onto 3.5 km s"' pixels. We add Gaussian noise to 
each spectrum^ assuming a flux-dependent root-mean-square (rms) 
amplitude given by"* (S/N)-' F(v) for S/N = 10, 30, and 50. Note 
that the typical S/N of STIS observations average S/N ~ 10 (e.g. 



3.1 Synthetic spectra 

Quite generally, the first step to compute a synthetic spectrum is to 
draw a random physical sightline from a simulation box at a given 
redshift z. A physical sightline is simply defined as the line between 



^ The redshift step is dz = 0.125. 

3 Unless stated otherwise, the S/N values quoted throughout this work de- 
note the S/J^ per pixel in the continuum. 
'' We assume a minimum, i.e. flux-independent noise level (S/N),„in = 10^. 
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Tripp et al. 2008), and hence we use our synthetic spectra with 
S/N = 10 when testing our results against observations in the next 
section. Spectra with S/N = 30 and S/N = 50 are used to assess 
the effect of the signal-to-noise ratio on the line statistics (see Ap- 
pendix A). Also, and even though a S/N of 50 is obviously much 
higher than the typical S/N of STIS observations, we choose it as 
a standard value for the spectra used in the analysis of the phys- 
ical conditions of O vi bearing gas in our simulation - presented 
in Sec. 5 - since this results in a larger sample of absorption lines. 
Besides, the Cosmic Origins Spectrograph (COS) recently installed 
on HST is expected to provide data at such high S/N for reasonable 
integration times^, although at somewhat lower spectral resolution. 

We fit O VI /1 103 1.93 A absorption features identified in our 
spectra with individual Voigt-profile components using a modified 
version of the package autovp (Dave et al. 1997), which includes 
the Voigt-profile approximation of Tepper-Garcia (2006). The fit 
results include the column density Nqvi and Doppler parameter 
foovi with corresponding uncertainties, as well as the rest equivalent 
width W r for each component. An example of a synthetic spectrum 
and its corresponding fit is shown in top panel of Fig. 1. 

The analysis of our 5000 spectra in the range < z < 0.5 
yields a total surveyed redshift path Az = 98.5, along which we 
identify a total of 1093, 2246, and 3115 Ovi absorbers for 
S/N=10, 30, and 50, respectively, thus obtaining a statistically sig- 
nificant sample of Ovi absorption lines. This allows a detailed 
analysis of the line parameter distributions for each absorption 
line sample obtained from spectra with different S/N, discussed 
in detail in Appendix A. Briefly, we show that we can reliably 
identify lines with Ovi column densities above log(Novi/cm"-) » 
12.3 at S/N = 50, log(Novi/cm-2) a 12.6 at S/N = 30, and 
logNovi/cm"^ a 13.2 at S/N = 10, even though we systematic- 
ally underestimate the total column density along each sightline, 
indicating that we miss a fraction of the absorption lines present in 
the spectra. For all identified lines, however, we are able to meas- 
ure rest equivalent widths quite accurately. Moreover, we find that 
the Doppler parameter distribution is sensitive to signal-to-noise 
ratios S/N < 10, but fitted Doppler parameters foovi 5 40km s"' 
are reliable if S/N > 30. Furthermore, we show in Appendix B 
that the line parameter distributions are converged with respect to 
box size and resolution. These results are of great importance when 
testing the predictions from simulations against observations us- 
ing quantities derived from synthetic spectra, since they allow us to 
bracket the sources of a potential disagreement between predicted 
and observed quantities more easily. Also, they are important for 
the analysis of the connection between physical conditions of the 
absorbing gas and direct observables, i.e., Ovi column density and 
Doppler parameter (see Sec. 5. 2). 

For ease of discussion, in what follows we shall refer to our 
Ovi absorber sample described above as sample 1. The particular 
line sample obtained from spectra with S/N = 10 is discussed in 
the next section where we compare the results from our simula- 
tion to observations. We generate and analyse another 5000 spectra 
with S/N = 50 for our fiducial run at z = 0.25. In these spectra 
we identify a total of 3034 O vi absorbers over a total redshift path 
Az = 97.7, and we shall refer to this O vi line sample as sample 
2. The latter is used and analysed in detail in Sec. 5. Note that the 
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Figure 1. The top panel shows a synthetic, normalised spectrum with 
S/N = 50 (black) and its corresponding fit (red) along a random sight- 
line through a simulation box of size L = 50/i"' Mpc at z = 0.25, showing 
a series of Ovi /11031.93 A absorption features. The vertical black marks 
flag the positions of the centroids of individual identified components. Note 
that only the relevant velocity range is shown. The middle and bottom pan- 
els show, respectively, the optical-depth weighted overdensity (magenta) 
and temperature (orange) along the sighthne in redshift-space (see Sec. 5). 
The values included in the middle panel coirespond to the optical-depth 
weighted gas temperature, optical-depth weighted overdensity, O vi column 
density, and O vi Doppler parameter for the line flagged by the blue dashed 
line. That the optical-depth weighted quantities correspond in each case to 
an optical-depth weighted average over the full line profile. Note in partic- 
ular the correspondence between density peaks and absorption features 



number of identified components and the total surveyed redshift 
path in sample 2 are very similar to the corresponding quantities of 
sample 1 . Also, the distribution of column density, Doppler para- 
meter, and rest equivalent width resulting from sample 2 are very 
similar to those of sample 1, but not identical, since the latter aver- 
ages the evolution of the O vi absorbers from z = 0.5 -^ 0, while 
the former is restricted to z = 0.25. 



Note, however, that our results are not sensitive to the assumed (S/N),^in 

value. 

^ So far, QSO spectra with a mean S/N = 50 have been obtained with COS 

for a total observing time of ~ 9 ks, compared to S/N = 10 for ~ 28 ks using 

STIS. 



4 COMPARISON TO OBSERVATIONS 

The Ovi line parameter distributions (presented in Appendix A) 
from our sample 1 for an adopted S/N = 10 can readily be 
compared to three well measured Ovi observables, namely, the 
column density distribution function (CDDF for short), the equival- 
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ent width distribution, and the correlation between Doppler para- 
meter and column density. 

The top-left panel of Fig. 2 displays the CDDF 
/(Novi) = '5^A/'/3Novifl/t' '•^■' the number of lines per logar- 
ithmic column density interval per unit absorption path length 
X (see Sec. 5.1.3). Comparison of the result from our fiducial 
run to observations by Thom & Chen (2008b, black symbols) 
shows that the amplitude of the predicted CDDF is slightly lower 
at column densities [10'^, 10'''] cm"^, but it agrees well with 
the high-column density data point, which has, however, a large 
associated uncertainty. Assuming the CDDF can be described 
in terms of a single power-law, Danforth & Shull (2008)* find 
that a fit in the column density range [10'', 6.3 x 10''*]cm"^ to 
the CDDF obtained from their Ovi absorber sample results in 
a power-law index /J = -1.98 ± 0.11. For the column density 
range [1.6x 10'^ 6.3x 10'^]cm-^ we findyS = -2.11 ± 0.15 
(solid line), which is consistent with the observed value. Note 
that this slope (with the appropriate normalisation) is also con- 
sistent with the data from Thom & Chen (2008b). If instead 
we consider only lines with column densities in the narrower 
range [1.6x 10'^ 2x 10'*]cm-^ we find /? = -1.78 ± 0.08 
(dashed line), which fits the distribution better as can be judged 
by the smaller uncertainty in p. We conclude that the slope 
of the Ovi CDDF resulting from our fiducial run is consistent 
with current observational constraints up to column densities 
Nqvi = 5 X 10'** cm"^, but the amplitude is lower (by a factor 2) 
than observed. It is worth noting at this point that a perfect match 
between predicted and observed CDDF would be surprising, since 
the stellar yields used in our simulation are uncertain at the factor 
of two level (see Wiersma et al. 2009b, their Appendix A3). 

Given that most of our identified lines lie on the linear part 
of the curve-of-growth (see Fig. A2), there is a (nearly) one-to- 
one correspondence between the predicted CDDF and the predicted 
equivalent width distribution, and we thus do not expect the latter to 
agree well with the observed one. The comparison between the cu- 
mulative equivalent width distribution both from observations and 
our sample 1 is shown in the top-right panel of Fig. 2. The ob- 
servations are taken from Danforth & Shull (2008, their Table 4, 
containing 75 confirmed equivalent width measurements) and from 
Tripp et al. (2008, their Table 2 with 77 single components). While 
Tripp et al. (2008) explicitly distinguish between components and 
systems, Danforth & Shull (2008) do not make a clear distinction 
between both types of absorbers. We follow the approach by Tripp 
et al. (2008) for a better comparison with their results. We define an 
absorption system as a simply connected region along the normal- 
ised, fitted spectrum with an O vi optical depth at each pixel above 
a given threshold t,i,. We adopt^ Tth = 0.097, which corresponds 
to the central optical depth of a Ovi -11031.93 A absorption line 
with a column density log(Novi/cm"^) = 12.3 - corresponding to 
the smallest detectable column density in our spectra with S/N = 
10; see. Fig. Al - and a Doppler parameter at the STIS resolution 
limit fo^in = 4.2 km s"'. Note that the resulting rest equivalent dis- 
tribution for systems is not completely insensitive to the adopted 
value of Tth. As can be seen in the top-right panel of Fig. 2, our 
fiducial run is not able to reproduce the amplitude of the observed 

* We note at this point that several studies (e.g., Wakker & Savage 2009; 
Cooksey et al. 2010; Williger et al. 2010) have disputed a subset of the line 
identifications of Danforth et al. (2006) and Danforth & Shull (2008), of 
which we make use in this section for comparison to our results. 
^ For comparison, Tripp et al. (2008) estimate a detection threshold given 
by Tth = 0.1 for their data with an average S IN x 10. 



distribution, neither for components nor for systems, as anticipated 
given the disagreement between the amplitude of the observed and 
predicted CDDF. Also, our predicted turn-over at W^ ~ 100 mA in 
the distribution of components is stronger than observed. 

To gain a better insight into the reason behind the discrep- 
ancy between the predicted and observed equivalent width distri- 
butions for single components, we re-plot the predicted and ob- 
served distributions shown in the top-right panel of Fig. 2 in dif- 
ferential form, as shown in the bottom-left panel of the same fig- 
ure. Note that the data and the results from our spectra have been 
binned using A log W^ = 0.2 dex, and that the resulting distribu- 
tions have been normalised to unit area. It is apparent that the ob- 
served distributions are broadly consistent with each other, even 
though the Tripp et al. (2008) distribution peaks at higher W^. 
In contrast, the distribution resulting from our fiducial run shows 
an excess of absorption lines at W,- ~ 30 mA, as well as a lack 
of absorption lines with Wr > 100 mA, which explains why our 
predicted turn-over at W,. ~ 100 mA in the cumulative equivalent 
width distribution is stronger than observed. As is shown in Ap- 
pendix A these mismatches are not due to fitting inaccuracies but 
are rather intrinsic to our simulation. Note further that, while Tripp 
et al. (2008)'s measurements are restricted to W^ > 30 mA, Dan- 
forth & Shull (2008) estimate their O vi system sample to be com- 
plete down to W r ~ 10 mA. 

In principle, there are two plausible reasons for the lack of 
strong (i.e., higher equivalent width) lines. On the one hand, it 
could be due to the apparent deficit of absorption components with 
Novi > lO'^'^cm"^ in the predicted CDDF (top-left panel of Fig. 2), 
when compared to the extrapolated power-law, assuming of course 
that the latter gives a correct description of the data. On the other 
hand, this lack could also arise due to the absence of lines with 
Novi ~ lO'^'^cm"^ and Doppler parameters much larger than the 
predicted median fc-value(i'o VI = 12.9 kms"' at S/N = 50). Indeed, 
since an O vi line saturates at these high column densities, the equi- 
valent width is particular sensitive to the Doppler parameter. An 
Ovi line with Nqvi ~ 10'"'-'cm"^ can, for example, have an equi- 
valent width as high as W^ ~ 400mA if feovi ~ 50 kms"'. Lines 
with such high fc-values are, however, almost absent from our spec- 
tra (see Fig. A2). Since there is no reason to assume that the ob- 
served power-law can be extrapolated to high column density val- 
ues, and given the good agreement between our predicted and the 
observed CDDF slope, we consider the absence of lines with large 
fe-values at a given column density as the most plausible explana- 
tion for the lack of high equivalent width lines in our simulation as 
compared to the data. 

To further investigate this possibility, let us now consider the 
correlation between Doppler parameters and Ovi column dens- 
ities. In the framework of a unifying model, Heckman et al. 
(2002) have shown that a iovi - Nqvi correlation, as is observed 
in Ovi absorption systems in a variety of environments (Milky 
Way disk and halo. High- Velocity Clouds, the Magellanic Clouds, 
IGM), naturally arises in diffuse, radiative cooling gas at temperat- 
ures T ~ 10^ - 10* K. Nevertheless, there is still some controversy 
about the validity of such a model for IGM O vi absorbers. Indeed, 
while Danforth et al. (2006) find no such correlation in their data, 
Lehner et al. (2006) show that their (own and compiled) data fol- 
low the predicted relation rather well. Similarly, Tripp et al. (2008) 
find a correlation between the line widths and column densities in 
their O vi sample, but these authors note that the significance is not 
enough to support the model by Heckman et al. (2002). 

Bearing this controversy in mind, we compare the fcovi - Nqvi 
correlation resulting from our simulation to observations. For this 
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Figure 2. Comparison between the observed O vi line parameter distributions and results from our fiducial run using spectra with S/N = 10 (our sample 1), 
which roughly corresponds to the average signal-to-noise of the data. Top left: O vi column density distribution function. The A:-bars indicate the bin size, while 
the r/-error bars show Poisson single-sided lo" confidence limits based on the tables by Gehrels (1986). The dashed and solid lines show, respectively, a fit in 
form of a power-law to the CDDF in the column density ranges [1.6 X 10 '"'j 6.3 X 10 '''] cm"^ and [1.6 X 10 '^, 2x10 ^'^} cm"^. Top right: Cumulative line 
number density from our simulation for single components (red) and systems (blue), and observations (symbols). In each case, the black solid line shows the 
actual value at each given limiting equivalent width, while the blue and red lines display Poisson noise. Bottom left: Dift'erential column density distribution 
for single components. The result for our simulation is highlighted by the shaded area. Bottom right: Line Doppler width v.s. O vi column density from our 
synthetic spectra (red) and observations by Tripp et al. (2008, black) and Danforth & Shull (2008, green). For ease of comparison, our predictions and the 
data have been binned in column density; in each bin, the dots indicate the median value, and the x- and y-bars indicate the bin size, and the 25th and 75th 
percentiles, respectively (see text for details). 



purpose, we use the results by Danforth et al. (2006, their 40 O vi 
single component sample) and Tripp et al. (2008, their high quahty 
sample of 77 intervening O vi absorbing components) only, since 
both these samples are large enough, since each of them provides 
evidence either against or in favour of the predicted correlation, and 
because we use one of these data sets for the comparison between 
our predicted and the observed cumulative equivalent width distri- 
bution. To facilitate the comparison between our results and obser- 
vations, we bin the O vi Doppler parameters from our simulation 
and the data in column density bins of size AlogNovi = 0.3 dex, 
and compute the median, and the 25th and 75th percentiles in each 
bin. The result is shown in the bottom-right panel of Fig. 2. Appar- 
ently, the predicted correlation between Doppler parameters and 
column densities matches the observations at the lowest column 
densities log(Novi/cm"") < 13.5, but does not so at higher column 
densities. While the observations show an increase in the Doppler 
width with increasing column density, the result from our simula- 
tion shows no clear trend. We warn, however, that the trend sug- 



gested by the binned data is to be taken with caution, given the 
enormous scatter in the observations. 



The bottom-right panel of Fig. 2 clearly shows that 
the observed fo-values are much larger at column densities 
Novi > 5 X 10'^ than the corresponding median fc-values of our 
identified O vi absorbers. This is consistent with our suspicion that 
the mismatch between the predicted and observed equivalent width 
distribution is due to a lack of O vi lines with high column densities 
and large Doppler parameters. Furthermore, this result also indic- 
ates that the observed O vi absorption systems are subject to a sub- 
stantial broadening mechanism, either in the form of (small-scale) 
turbulence (which is not captured/properly modelled by the SPH 
scheme in our simulation) or Hubble broadening. The alternative 
that the observed Ovi absorbers arise in gas at higher temperat- 
ures, thus resulting in lines with larger Doppler parameters, can be 
ruled out, since a Doppler parameter of, e.g. , io vi = 40 km s"' cor- 
responds to a gas temperature T ~ 1.5 x 10* K for oxygen, which 
would result in a too low O vi ion fraction (novi/"o) < 10"", either 
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in CIE (Sutherland & Dopita 1993) or non-equilibrium conditions 
(Gnat & Sternberg 2007). 

It is noteworthy that a similar disagreement between the pre- 
dicted and the observed O vi cumulative line-number density and 
bo VI - No VI correlation is reported by Oppenheimer Sl Dave (2009) . 
These authors argue that turbulence is the crucial mechanism lead- 
ing to the large fe- values measured in IGM O vi systems, and solve 
the discrepancy between their simulations and observations by 
adding sub-resolution turbulent broadening as a function of hy- 
drogen density ~ partly constrained by observations- to the Dop- 
pler parameter of their simulated absorption lines. Quite remark- 
ably, they find that this approach, besides reproducing the b - Nqvi 
correlation by construction, simultaneously brings their predicted 
equivalent width distribution into better agreement with observa- 
tions. Since the CDDF slope resulting from our simulation is con- 
sistent with the data, while the fc-value distribution is not, it is quite 
plausible that broadening our O vi lines following the approach by 
Oppenheimer & Dave (2009) would help to reconcile our results 
with observations. Although this is tempting, we will not consider 
it here. Instead, we will present and discuss a series of modifica- 
tions to our reference run {i.e. simulation runs with different para- 
meters) and post-run variations thereof - including the addition 
of sub-resolution turbulence- and their implications for Ovi ab- 
sorbers statistics in a future paper. For now, we will proceed with 
the analysis of the physical properties of the O vi bearing gas in our 
reference run as is. 



5 PHYSICAL PROPERTIES OF O vi ABSORBERS 

In this section we analyse the physical properties (density, temper- 
ature, metallicity, ionisation state, baryon content) of the O vi bear- 
ing gas traced by the absorption features identified in our synthetic 
spectra. For the sake of simplicity, we restrict ourselves to the ana- 
lysis of our fiducial run at z = 0.25, which roughly corresponds 
to the median redshift of the O vi observations found in the liter- 
ature. Hence, in what follows, we will use our O vi line sample 2 
introduced at the end of Sec. 3.1. 



5.1 Optical depth-weighted physical quantities 

Defining physical properties such as the density or temperature of 
the gas responsible for the absorption features identified in simu- 
lated spectra is not a trivial task, and different methods have been 
described in the literature. Oppenheimer & Dave (2009, their equa- 
tion 7), for example, ascribe physical properties to simulated Hi 
and O VI absorption systems by weighting the desired quantity by 
the product of the SPH particle mass, the mass fraction of the cor- 
responding element {e.g. oxygen), and the ionisation fraction of 
the corresponding species {e.g. Ovi) at the line centre, taking both 
peculiar and thermal velocities into account. Schaye et al. (1999), 
on the other hand, define the density (temperature) of the absorb- 
ing gas at the line centre as the sum of the density (temperature) 
of all gas elements that contribute to the absorption in that pixel, 
weighted by their contribution to that pixel's optical depth. A sim- 
ilar approach is followed by Richter et al. (2006) to ascribe a dens- 
ity and a temperature to the absorption features identified as Broad 
Lya Absorbers (BLAs) in their simulations. 

By definition, optical-depth weighted physical quantities dir- 
ectly relate an absorption feature identified in a synthetic spectrum 
to the absorbing gas in the simulation, and this gives insight about 
the physical state of the gas in which absorption actually takes 



place. This in turn allows for a meaningful comparison between 
quantities obtained from a spectral analysis {e.g. , column density, 
Doppler parameter) and the 'true' physical properties as given in 
the simulations {e.g., volume density, temperature). For example, 
for a uniform density, uniform temperature gas cloud, the optical- 
depth weighted properties of the cloud would be simply its given 
density and temperature. In case there is, e.g., a density (temper- 
ature) gradient towards the centre, then naturally the denser gas 
will be more important in determining the absorption line prop- 
erties. The optical-depth weighted quantity properly reflects this, 
and provides the properly weighted average density or temperat- 
ure across the absorbing cloud. We therefore choose to use optical- 
depth weighted physical quantities for our analysis of simulated 
spectra. In contrast to the approaches mentioned above, however, 
we do not only consider the line centre, but compute a weighted 
average over the line profile. So, to compute a physical quantity, 
say the density, associated with a certain line, we first compute the 
optical-depth weighted density in redshift space along the sightline 
as in Schaye et al. (1999). Next we compute the average of the 
optical-depth weighted density over the line profile, weighted by 
the optical depth of each pixel and assign this last weighted aver- 
age to the line. We have compared the physical quantities resulting 
from this approach with using just the value of the optical-depth 
weighted quantity at the line centre and found that the difference is 
small, with averaging over the full profile typically giving slightly 
smaller values. 

In what follows, quantities weighted by Ovi optical depth 
will be denoted by adding a corresponding subscript; thus, for ex- 
ample, the optical-depth weighted gas overdensity (temperature) is 
denoted by Aqvi (7'ovi)- Note that we define A = Ph/ (Pb), where 
(Pb) = 4.18 X 10"^' {h/0.13)^ g cm"^ is the mean cosmic baryon 
density. An example of the optical-depth weighted temperature and 
overdensity along a random sightline through a simulation box at 
z = 0.25, as well as the corresponding quantities for a given Ovi 
absorption line are shown in middle and bottom panels of Fig. 1 . 



5.1.1 Temperatures and overdensities 

In a recent study, Wiersma et al. (2009b) investigated the distri- 
bution of metals using a simulation run which assumed the same 
reference model as our fiducial run but at 8x lower mass resolu- 
tion and twice as large a box size. Their results, which were shown 
to be insensitive to box size and nearly converged for both their 
and our fiducial runs, show that the diffuse, photo-ionised IGM, 
i.e. , gas at overdensities A < 10'-' and temperatures 10' K < 
T < 10"*^ K, harbours a large fraction of the mass in the simu- 
lation, but contains a negligible fraction of the metals. They find 
that the metals are mostly spread over low-density structures at 
temperatures T > 10^ K, despite the inclusion of metal-line cool- 
ing, indicating that the Warm-Hot Intergalactic Medium (WHIM) 
at 10^ K < r < 10'' K and A < 10' contains a significant amount 
of metals (Wiersma et al. 2009b, their Figure 10). As discussed by 
these authors, the reason that a very high fraction of the metals are 
in low-density gas at T > 10^ K is that high- velocity winds trans- 
port metals from galaxies out to large distances, and these winds 
shock-heat the gas to such high temperatures. 

Motivated by their findings, here we want to investigate if and 
how the baryon and metal distributions, as well as different gas 
phases in our simulation are traced by Ovi detected in absorp- 
tion. To this end, we compute an optical-depth weighted temper- 
ature, To VI 5 and overdensity, Aovi> for each of our Ovi absorbers 
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in sample 2, and compare the resulting (rovi, Aovi)-distribution to 
the intrinsic mass and metal distributions in our simulation. 

The top panel of Fig. 3 shows the gas mass distribution 
"J^^gas/C^gas.totSlog AiSlogT) (coloured areas) at z = 0.25 
and the distribution of Ovi absorbers (white contours) in the 
temperature-density plane. For the case of the absorbers (contours) 
the axes correspond to the O vi optical-depth weighted temperature 
and density. The colour coding shows the amplitude of the gas mass 
distribution. The red/green/blue/black areas enclose 50/75/90/100 
per cent of the total gas mass, with the gas mass fraction in each 
bin given by the value indicated in the colour bar. The white con- 
tours contain, starting from the innermost, 25, 50, 75, 90, and 99 per 
cent of the total number (3034) of identified O vi absorbers. For ref- 
erence and for the subsequent analysis, we include the sub-panels 
below and to the right which show, respectively, the distributions 
(normalised to their peak value) marginalised over A (lower sub- 
panel) and T (right sub-panel), of gas mass (red), and O vi absorbers 
(black). Note that we have excluded from this and all subsequent 
phase diagrams the star-forming gas (i.e., the interstellar medium) - 
defined as gas with densities exceeding our adopted star-formation 
threshold n^ = 0. 1 cm"^ which corresponds to* A « 3 x 10' at 
z = 0.25 - since the temperature of this gas simply reflects the pres- 
sure, imposed via an equation of state P oc p"*'', of the unresolved 
multiphase ISM (see Schaye & Dalla Vecchia 2008, for details). 

Focusing first on the gas mass distribution, we can see that 
approximately half of the baryonic mass in our simulation is 
contained in a phase corresponding to the diffuse IGM (red re- 
gion), which shows a tight correlation between temperature and 
(over)density governed by the balance between adiabatic cooling 
due to the cosmic expansion and photo-heating by the meta-galactic 
UV background (Hui & Gnedin 1997). We also find a signific- 
ant amount of gas at low overdensities (1 < A < 10") and much 
higher temperatures (10' K < T < 10' K), corresponding to shock- 
heated material. The plume of the distribution at A > 10' and 
T > 10' K , represents hot gas in galaxy clusters, i.e. the so-called 
intra-cluster medium (ICM). Note that there is no well-defined de- 
marcation between gas phases at different temperatures and over- 
densities (see red histograms); rather, the transition between these 
phases is smooth. 

In contrast to the extended temperature and overdensity distri- 
butions of the gas mass in our simulation, the distribution of Tovi 
and Aovi shown by the white contours is constrained to a relat- 
ively narrow temperature range 10*' K < T < 10' K and overdens- 
ities 1 <K A < 10' . Our simulation thus indicate that Ovi traces 
shock-heated material with temperatures around T ~ lo"*"' K, 
and at slightly higher overdensities than the diffuse IGM, around 
A ~ 10 - 10^ typically. The comparison to the overall gas mass 
distribution also suggests that Ovi traces a significant amount of 
the baryons in the simulation. We will discuss the baryon content 
of O VI bearing gas in more detail in Sec. 5.1.3, and will show that, 
contrary to our expectation, the Ovi bearing gas contains only a 
small fraction of the warm-hot baryons. 

The middle panel of Fig. 3 shows the gas metal mass distri- 
bution d~ Mz/(Mz,tot d log A d log T). Each particle's metal mass is 
defined by niz = Z^n, x m^, where Zj^, is the smoothed metallicity 



The relation between hydrogen number density nn and (baryonic) over- 
density A is given by 
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Figure 3. Distribution (coloured areas) of gas mass (top), gas metal mass 
(middle), and Ovi mass (bottom) in the temperature-overdensity plane at 
z = 0.25, together with the distribution of O vi absorbers (white contours) of 
sample 2. The histograms at the bottom and to the right of each panel show 
the corresponding gas (red) and O vi absorbers (black) distributions mar- 
ginalised over T and A, respectively, normalised to their peak values. The 
colour coding has been chosen in such a way that the red/green/blue/black 
areas contain 50/75/90/100 per cent of the total gas mass (top), total gas 
metal mass (middle), total Ovi mass (bottom), respectively. The white con- 
tours enclose, starting from the innermost, 25, 50. 75, 90, and 99 per cent 
of the total number of O vi absorbers. Note that the white contours and the 
black histograms are the same in all three panels. 
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(see Sec. 5.1.2 and Wiersma et al. 2009b, for a detailed discussion 
on smoothed metallicities; we note that this plot would look nearly 
identical if particle metallicities had been used). Again, the col- 
our scheme shows the amplitude of the distribution with the colour 
cuts chosen in such a way that the red/green/blue/black areas en- 
close 50/75/90/100 per cent of the total gas metal mass, with the 
gas metal mass fraction in each bin given by the value indicated in 
the colour bar. The white contours are the same as in the top panel. 
Clearly, we can distinguish two phases that harbour a substantial 
fraction of the metals in the simulation at z = 0.25: shock heated 
material at temperatures T > 10' K and overdensities A < 10^, 
and cold-warm gas at T < 10"*' K and A > 10^. In agreement with 
Wiersma et al. (2009b, their Fig. 7) we estimate that these phases 
together contain ~ 30 per cent of the total metals in the simula- 
tion, while the rest is contained in star- forming gas (~ 10 per cent) 
or locked up in stars (~ 60 per cent). We find that nearly 20 per 
cent of the metals in our simulations are contained in WHIM gas 
at 10' K < r ^ 10' K, but the O vi absorbers together contain only 
< 0.1 per cent of the metal budget. Hence, O vi absorbers do not 
trace the bulk of metals, even though the majority arise in gas with 
a high metal fraction, as can be judged from the (rovi, Aqvi) dis- 
tribution. 

Presumably, there is a continuous exchange of material 
between the cold-warm, diffuse IGM and these metal-rich gas 
phases. Intergalactic gas is accreted onto the potential wells of 
galaxies and thereby shock-heated. At the same time, galactic 
winds shock-heat and transport metal-rich gas from the high- 
density, star-forming regions into the low-density IGM, some of 
which may then cool down and may eventually be re-accreted onto 
galaxies to fuel star formation ("intergalactic fountain"). In this 
scenario, and given the location of the O vi distribution with respect 
to the gas metal mass distribution (middle panel), Ovi absorbers 
would mainly arise in enriched, shock-heated galactic-wind mater- 
ial that has probably started to mix with the ambient gas and cool 
down, wandering from the high-temperature, low-density region to 
the low-temperature, high-density regime on phase space. As such, 
many of the Ovi absorbers tracing warm-hot gas may be short- 
lived. We will provide evidence for this later on when we discuss 
the metallicity of Ovi bearing gas below and in Sec. 5.1.2. 

To better understand the nature of these O vi absorbers, we 
investigate if (and how much of) the O vi bearing gas in our sim- 
ulation is traced by the Ovi we see in absorption. To this end, 
we compare in the bottom panel of Fig. 3 the O vi mass distribu- 
tion d^ Mom/(Movi,tot d log A d log T) to the distribution of O vi ab- 
sorbers in the (Tovi, Aqvi) plane. A particle's oxygen mass is given 
by mo VI = (tiovi/no) x Xq x mg, where Xq is the smoothed oxygen 
mass fraction. Colours and contours have the same meaning as in 
the middle and top panels. We find that the bulk of the O vi is found 
in gas at moderate overdensities A ~ 10"' - 10^', and temperat- 
ures between T = lO'^K and 10^ K, which is consistent with the 
overall distribution of our O vi absorbers. Also, the marginalised 
temperature distribution (lower sub-panel, red histogram) is clearly 
bi-modal, showing that O vi is distributed among two phases with 
different temperatures: T ~ lO'*' K, corresponding to photo-ionised 
gas, and T ~ 10" K, corresponding to coUisionally ionised gas. We 
find that more than 65 per cent of the O vi seen in absorption traces 
mainly the hot gas phase at temperatures T > 10' K. This means 
that approximately one third of the O vi absorbers arise in cooler, 
photo-ionised gas. This number agrees with the result of Tripp et al. 
(2008) who find compelling evidence that 34 per cent of the O vi 
absorbers in their sample arises in gas at T < 10' K. 

Our finding that most of the Ovi arises in hot gas is not 



caused by the limited S/N of our spectra, nor by the limited res- 
olution, since qualitatively the same result is found in noise-free 
spectra with very high resolution. As can be judged by the amp- 
litude of the (two-dimensional) distribution around T ~ 10**' K 
and 10 < A < 10^, the Ovi content in this gas phase is low, which 
could mean that also the O vi column density is too low for all of 
this gas to be detectable in absorption. Our assumption is reinforced 
by the fact, going back to the middle panel, the overall metal con- 
tent of gas at temperatures and overdensities comparable to the low 
temperature Ovi bearing gas is quite low. The fact that the peak 
temperature of the hotter phase is slightly shifted to lower values as 
compared to the high-temperature peak of the intrinsic temperat- 
ure distribution, suggests that O vi absorption arising in cooler gas 
phases overlaps in redshift space with Ovi absorption from hot- 
ter gas, thus leading to slightly lower optical-depth weighted gas 
temperatures. To a lesser extent, it is a consequence of averaging 
optical-depth weighted quantities over the full line profile, as dis- 
cussed in Sec. 5.1. 

An interesting fact is that the low temperature O vi bearing gas 
closes the gap between the high-temperature, low-density and the 
cold, high-density metal-rich gas phases shown in the middle panel, 
which suggests that some of the gas traced by Ovi has already 
cooled down to an equilibrium temperature T ~ 10**' K, which 
is typical for enriched gas at z = 0.25 (see Sec. 5.3, Fig. 10). 

Summarising, we find that, while the mass in our simulation is 
more or less equally distributed between two phases with temperat- 
ures T < 10"' K and r > 10"' K, and overdensities 0.1 < A < 10^ 
Ovi traces gas at high temperatures T > 10' K and moderate over- 
densities 10 < A < 10^, which are typical for the warm-hot inter- 
galactic medium (WHIM) in our (Wiersma et al. 2009b) and other 
simulations (Cen & Ostriker 1999; Dave et al. 2001; Bertone et al. 
2008; Tornatore et al. 2009). This gas phase typically contains ~ 40 
of the baryons in our simulation, but the O vi bearing gas turns out 
to contain only a small fraction of the cosmic baryons at low red- 
shift (see Sec. 5.1.3). Furthermore, we find that roughly 40 per cent 
of the metals in our simulation are distributed between warm-hot 
gas at moderate to high overdensities 1 < A < 10' (~ 30 per 
cent), and cold star-forming gas at very high densities (~ 10 per 
cent). While Ovi certainly arises in gas containing metals, it is not 
a tracer of the main metal reservoirs in our simulation. 

Our result that absorption by O vi in our simulations is pref- 
erentially related to gas at temperatures 10**' K < T < 10* K, with 
a high fraction (more than 65 per cent) having temperatures fall- 
ing in the range 10' K < T < 10* K, i.e. , in the low temperature 
regime of the WHIM, is quite interesting, since Oppenheimer & 
Dave (e.g. 2009) find exactly the opposite trend, namely that the 
overwhelmingly majority of the O vi in their simulations is found 
in gas a temperatures T ~ lO'*-^*''-^ K, typical of photo-ionised gas. 
It is important to note that our simulation does show a gas phase 
containing O vi at photoionisation temperatures (see Fig. 3, bottom 
panel), which is however not detected in absorption, most probably 
due to its low overall metal content. In Sec. 5. 3 we will present a 
deeper analysis of the disagreement between the findings by Op- 
penheimer & Dave (2009) and our results, which we will attribute 
to their neglect of the effect of photo-ionisation on metal-line cool- 
ing. 



5.1.2 Metallicity 

In this section, we analyse the local metallicity of the gas traced by 
the O VI absorbers (in sample 2). We proceed in the same manner 
as with gas temperatures and overdensities, and estimate optical- 
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Figure 4. Top row: Gas mass distribution (coloured areas) on the overdensity-metallicity (left) and temperature-metallicity (right) planes, together with the 
distribution of the O vi absorbers (white contours) for our fiducial run REF_L050N512 at z = 0.25. Note that the marginalised overdensity distributions (lower 
sub-panels) are identical to those in the top panel of Fig. 3 (right sub-panel). Bottom row: Gas metal mass distribution on the overdensity-metallicity (left) 
and temperature-metallicity (right) planes for our fiducial run REF_L050N512 al z = 0.25. White contours are as in the top row. Note that the marginalised 
temperature distributions (lower sub-panels) are identical to those in the middle panel of Fig. 3 (lower sub-panel). Colour coding and contour levels are as in 
Fig. 3. We use a solar metallicity of Zq = 0.0127 to express the absolute metallicities predicted by the simulation in solar units. 



depth weighted metallicities Zqvi for each of our identified O vi ab- 
sorbers. Note that an accurate measurement of the elemental abund- 
ances in the absorbing gas is very important, since this quantity is 
crucial for an indirect estimate of the baryon fraction in the gas 
traced by a given absorption species (see Sec. 5.1.3). 

In our simulation, as is done in all OWLS runs, we trace how 
stellar evolution alters the abundance of 1 1 elements (H, He, C, N, 
O, Ne, Mg, Si, Fe, Ca, S) explicitly, while simultaneously follow- 
ing the metallicity of each particle (see Wiersma et al. 2009b, for 
further details). We choose to use smoothed metallicities - as op- 
posed to particle metallicities - for the various processes that are 
metallicity-dependent, in particular for the calculation of radiative 
metal cooling. The smoothed metallicity of a particle is defined 
as the ratio of the SPH estimates of the metal mass density and 
the total gas mass density at the location of the particle, while the 
particle metallicity is given by the ratio of the particle's metal mass 
and its total gas mass. Wiersma et al. (2009b) have shown that 
in low-metallicity gas, smoothed metallicities are generally higher 
than particle metallicities. The use of smoothed metallicities in part 
counters (but does not solve) the lack of metal mixing inherent 
to SPH, which leads to differences when compared to the use of 
particle metallicities (although these differences decrease with in- 
creasing resolution). First, the use of smoothed metallicities results 



in higher fractions of metals residing in gas both at lower temperat- 
ures (r < 10' K as compared toT ~ 10* K) and lower metallicities. 
Second, the use of smoothed abundances enhances radiative cool- 
ing and increases the predicted SFR. These differences may be im- 
portant for the comparison of our results to the results from studies 
employing particle metallicities (e.g. Oppenheimer & Dave 2009). 
In the top row of Fig. 4, we show the gas mass distribu- 
tions S^Mjas/CMuas.totSlogZcJlog A) (left panel; coloured areas) 
and (S^Mjas/CMjas.tot'JlogZcJlogr) (right panel, coloured areas), 
together with the distribution of the O vi absorbers in the (A, Z) 
and (T, Z) planes (white contours), respectively. Again, for the 
case of the absorbers the axes correspond to the O vi optical-depth 
weighted temperature Tovi, density Aqvi, and metallicity Zovi. The 
meaning of colours and contours is the same as in Fig. 3. The his- 
tograms in the lower and right sub-panels show the correspond- 
ing distributions marginalised over metallicity and overdensity (left 
panel) or temperature (right panel) normalised to their peak values. 
The panels in the bottom row show the gas metal mass distribu- 
tions d^ MzKMz.tot d log Z d log A) (left panel; coloured areas), and 
d^ Mz/(Mz,ot 51ogZ31og T) (right panel; coloured areas), as well 
as their corresponding marginalised distributions (histograms in the 
low and right sub-panels), together with the distribution of the O vi 
absorbers (white contours). Note that these panels follow the same 
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colour scheme and include identical contours as the corresponding 
top-row panels, but that the !/-axis range is smaller. 

Several interesting conclusions can be drawn from this figure. 
First, the top-row panels reveal that a high fraction of the baryonic 
mass in our simulation is contained in gas with local metallicit- 
ies -2 < log(Z/ZQ) < -1. In contrast, Ovi absorbers clearly trace 
over-enriched material with metallicities log(Zovi/Zo) > -1. The 
medianmetallicity of thegas tracedby OviisZovi ~ 0.6 Zq, which 
is almost an order of magnitude higher than the median gas metal- 
licity at similar overdensities and temperatures. Note that the use 
of optical-depth weighted metallicities ensures that we track (and 
average over) the metallicity of all individual particles that contrib- 
ute to the Ovi absorption. Hence, Zqvi is an estimate of the local, 
i.e. around the spatial resolution limit, metallicity. Note that the 
overall mean metallicity of the structures giving rise to the O vi ab- 
sorption - as is derived from observations by using the ratio of the 
O VI column density to the total H i column density in the absorbing 
gas - can be much smaller than Zqvi- In our case, this is supported 
by the fact that gas at overdensities and temperatures comparable 
to the centres of mass of the (Aqvi, Zqvi) and (rovi, -Zovi) distri- 
butions display a large range of metallicities. The large scatter in 
metallicities, in particular at the lowest overdensities, thus suggests 
that the metal distribution is quite inhomogeneous. 

In order to investigate more closely the metal distribution in 
the structures where the O vi absorption takes place, we proceed as 
follows: First, we fit all Hi /11215.67 A absorption features along 
the same spectra from which we obtain our O vi sample 1 . Next, we 
identify every O vi component in the latter sample that is coincid- 
ent with a given H i component within ± 10 km s" ' in velocity space. 
In this way we obtain a sample of well-aligned H i - O vi absorber 
pairs. Assuming that the coincidence in velocity space implies that 
the H I and O vi absorption arise in the same gas, we estimate the 
(average) metallicity Z(0 vi) of the structure giving rise to the ab- 
sorption via 
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We compute Z(Ovi) for each well-aligned absorber pair us- 
ing the corresponding optical-depth weighted ionisation fractions, 
('Jhi/«h)hi and («ovi/no)ovi, and the Ovi optical-depth weighted 
hydrogen fraction, (Xh)ovi- Also, we compute for each individual 
H I- and O vi-member of a well-aligned absorber pair a correspond- 
ing pair of optical-depth weighted metallicities, Zhi and Zqvi, re- 
spectively. The distributions of the resulting metallicities for all 
well-aligned H i and O vi components in our sample are shown in 
Fig. 5. Clearly, the distribution of the mean metallicity estimated 
from the Ovi -Hi column density ratio, Z(Ovi), shows a larger 
scatter and has a peak value which is lower by an order of mag- 
nitude than that of the Zqvi distribution. In contrast, the Z(Ovi) 
distribution is comparable in width and range to the distribution of 
Zhi, although it peaks at a slightly lower value. Interestingly, there 
seems to be a sub-population of well-aligned Hi- Ovi absorber 
pairs with metallicities around 0.6 Zq, where the optical-depth 
weighted metallicity of the H i-member matches well the optical- 
depth weighted metallicity of the O vi-member. 

This exercise readily shows that local gas metal content, as 
measured by the Ovi optical-depth weighted metallicity, can dif- 
fer substantially from the overall, average metallicity of the struc- 
ture giving rise to the observed absorption. This in turn strongly 
suggests that the metallicity distribution is highly in-homogeneous, 
and that O vi absorbers are biased towards high-metallicity patches 
embedded in gas with an otherwise lower, average metallicity. 



Ovi weighted 
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Figure 5. Metallicity distribution for our sample of well-aligned H i - O vi 
absorber pairs. The red (blue) histogram shows the distribution of O vi (H i) 
optical-depth weighted metallicities. The distribution of metallicities estim- 
ated from the the Ovi - Hi column density ratio (Nqvi/Nhi) is shown by 
the black curve (see eq. 1 and corresponding text for details). 



We note that the patchy metal distribution may in part be due 
to the metal-mixing problem inherent to SPH which is alleviated 
(but by no means solved) by the use of smoothed abundances, 
as mentioned above. While the implementation of metal diffusion 
(Greif et al. 2009; Shen et al. 2009) could, depending on the choice 
of parameter values, reduce the scatter in the metallicities, this may 
in fact be in conflict with observations which indicate that inter- 
galactic metals are poorly mixed on small scales (Schaye et al. 
2007). Our simulations, on the other hand, seem to agree with the 
latter observational fact, at least in terms of well-aligned H i - O vi 
absorber pairs. We will present a detailed analysis of the alignment 
statistics, as well as of the physical properties of aligned pairs, in a 
companion paper {Tepper-Garcia et al, in preparation). 

In terms of the metal mass distribution, we can see in Fig. 4 
that the vast majority of the metals in our simulation reside in gas 
with metallicities -1 < log(Z/Zo) < 0, temperatures T > lO'' K, 
and overdensities 10"' < A < 10^. From the overlap of the Ovi 
distribution with the overall metal mass distribution, we may con- 
clude that the O vi we see in absorption traces the moderate-density, 
warm, enriched gas relatively well. Nevertheless, despite the coex- 
istence of the O VI absorbers with this gaseous phase on the (T, Z)- 
and (A, Z)-planes, the gas traced by O vi contains actually only a 
vanishingly small (less than 0. 1 per cent) fraction of the metals (see 
Sec 5.1.1). 

To sum up, the results of this section indicate that O vi traces 
over-enriched gas, which is neither representative for the bulk of 
the baryons nor for the bulk of metals. The metallicity of the gas 
traced by O vi is high, with a narrow distribution extending over the 
range (0.1, 1) Zq, and a peak at Zqvi ~ 0.6 Zq, which is consistent 
with results from other simulations (e.g. Cen & Ostriker 2006; Op- 
penheimer & Dave 2009), and metallicity measurements in metal- 
line absorption systems at low redshift (e.g. Savage et al. 2002; 
Prochaska et al. 2004; Tripp et al. 2005; lenkins et al. 2005; Cook- 
sey et al. 2008). In spite of the general agreement with observations, 
we want to emphasise again that our quoted values refer to local 
estimates of the metallicity, while observed values are generally 
derived from metal-line systems aligned (i.e. within a given range 
in velocity space) with H i absorption, which can lead to signific- 
antly lower mean metallicities if the metal distribution is patchy 
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Figure 6. Distribution of optical-deptli weiglited Ovi fraction (top-left), 
local oxygen abundance (top-right), and local hydrogen mass fraction 
(bottom-left) for our sample of 3034 O vi absorbers at z = 0.25 in spec- 
tra with S/N = 50 (sample 2). The bottom-right panel shows the distri- 
bution of the baryon density, relative to the cosmic mean, traced by O vi 
along individual sightlines. Note that we use the default cloudy solar oxy- 
gen abundance (no/iH)o = 4-. 9 X lO"'' to compute [0/H]ovi- 



(see Schaye et al. 2007, for a discussion). We showed that this can 
be indeed the case by comparing optical-depth weighted metallicit- 
ies with metallicity estimates derived from O vi -H i column density 
ratios for well-aligned H i and O vi absorber pairs. 



5.1.3 Baryon content 

The baryon content of gas traced by O vi is a quantity of great in- 
terest since it has long been anticipated by cosmological simula- 
tions that nearly half of the baryons in the Universe (absent in ob- 
servational inventories, see e.g. Fukugita 2004) might be hidden in 
the WHIM (Cen & Ostriker 1999; Dave et al. 2001; Bertone et al. 
2008), and that this gas phase might be detectable through its dif- 
ferent emission and absorption signatures, in particular absorption 
by Ovi (Cen & Ostriker 1999). 

A standard approach used in the literature to indirectly estim- 
ate the baryon fraction of O vi bearing gas is to write 
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where ran is the hydrogen mass, pc the critical density, and Xh the 
(assumed) hydrogen mass fraction. The factor in brackets is the 
total surveyed physical length with Axi the so-called absorption 
path length (see e.g. Schaye 2001) along an individual sightline. 

One major caveat of this approach is that the estimated value 
for fib (Ovi) is highly sensitive to both the ion fraction and the 
metallicity of the gas and that, owing to the difficulties in measur- 
ing these quantities, plausible values need to be assumed. Usually, 
the Ovi ion fraction is set to (novi/"o) = 0.2 - corresponding to 
the peak ion fraction of Ovi at T ~ 2-3x 10' K either in CIE 
(Sutherland & Dopita 1993) or non-equilibrium (Gnat & Sternberg 
2007) - and the metallicity to [O/H] =0.1 dex (see e.g. Tripp et al. 
2000). 

Here we want to estimate the baryon content of the gas traced 



by our sample of O vi absorbers, and we want to do this in such a 
way that we can compare our result to the observed range of val- 
ues. Therefore, we follow the above approach, i.e., using equation 
(2), but exploit the advantage of our simulation, which allows us to 
measure the ion fraction (novi/«o)> oxygen abundance (no/iin), and 
hydrogen mass fraction Xh of the gas giving rise to each O vi ab- 
sorber. To this end, and to be consistent, we again use optical-depth 
weighted quantities, estimated for the O vi absorbers in sample 2. 
The resulting distributions of the values for each of these quantities 
are shown in Fig. 6 (top row; bottom-left). 

Inserting the values for each of the latter quantities and 
for each identified Ovi absorber in equation (2), we obtain 
fib(Ovi)/fii, = 0.014 over 5000 sightlines spanning a total ab- 
sorption path Ax = 139.3. The distribution of the baryon content 
along individual sightlines is shown in the bottom-right panel of 
Fig. 6. Using the same 5000 spectra', this time with S/N = 10 
rather than S/N = 50, we find fib(0 vi)/Qb = 0.0096. Either value 
is lower than values estimated from observations which range from 
fib(Ovi)/fib i 0.019 (Tripp et al. 2000), through ~ 0.045 (Thom 
& Chen 2008b), > 0.054 (Sembach et al. 2004), and 0.086 ± 0.008 
(Danforth & Shull 2008). 

Two things are noteworthy. First, as can be judged from the 
bottom-right panel of Fig. 6, the distribution of fib(Ovi) along in- 
dividual sightlines shows a large scatter, reflecting the effect of cos- 
mic variance. This indicates that measurements of the baryon con- 
tent using data taken along different sightlines might vary drastic- 
ally, even in case that the ionisation fraction and the metallicity 
are known to high accuracy, which they are usually not (although 
note that our sightlines span Az = 0.02 at z = 0.25). Second, and 
related to the previous point, is the fact that the median optical- 
depth weighted ionisation fraction resulting from our simulation 
is (novi/no)ovi ~ 0.052, i.e., is a factor four smaller than usu- 
ally assumed, and the median optical-depth weighted metallicity 
is [0/H]ovi ~ -0.19 dex or about six times larger than the 'canon- 
ical' 10 per cent solar-value. Also, the distribution of both these 
quantities shows a large dispersion (see top panels in Fig. 6). This 
in turn implies that baryon density measurements from observa- 
tions which assumed fixed values for the Ovi ionisation fraction 
and the metallicity may be off (most probably overestimated) by 
non-negligible factors, independent of the effect of cosmic vari- 
ance. This in turn might explain why we find an average baryon 
density which is lower than the various measured values. 

A second quantity we can easily compute and compare to 
measured values is the baryon content, relative to the critical dens- 
ity, of O VI in all identified absorbers. The expression used to com- 
pute this quantity is very similar to equation (2) and is given by 
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where otq is the oxygen atomic mass. This quantity is dir- 
ectly measurable, and has the advantage that it is independent 
of the metallicity and ionisation corrections. If we use the re- 
cwvered, i.e., fitted Ovi column densities, we get fio VI = 1.1 x lO"' 
for S/N = 50, and fiovi = 9.3 x 10"** for S/N = 10, down to 
Wr = 10 mA. If instead we use the true, total Ovi column density 
along each sightline we find fiovi = 1.2 x 10"'. Note that the latter 
value exactly matches (neglecting numerical round-off errors) the 
O VI mass density relative to the critical density obtained by adding 



' In this case we identify a total of 1044 O vi components rather than 3034. 
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up the O VI mass of each SPH particle in the simulation and dividing 
by the simulation volume. For comparison with our results, Thom 
& Chen (2008b) report fiovi = 1-7 ± 0.3 x 10"' for W, > 30 mA., 
while Danforth & Shull (2008) find novi = 4.9 ± 0.4 x 10"' for 
Wr > 10 mA. Both values are higher than ours, but this is expec- 
ted since the amplitude of our CDDF is also lower compared to 
observations. 



5.2 Physical conditions of the gas and O vi observables 

Now that we have analysed the physical conditions (density, tem- 
perature, etc.) of the gas giving rise to the O vi absorption in our 
simulation, we want to make a connection between these and the 
actual O VI observables, i.e. the Ovi column density, Novi> and the 
O VI line width as measured by the Doppler parameter, Z?ovi- 

For this purpose we bin the set of (fcovi, Novi) values for our 
Ovi sample 2 using Aiovi = 2.0 kms"' and AlogNovi = 0.2 dex, 
and compute in each cell the median of the physical quantity un- 
der consideration, e.g. density. The result of this exercise is shown 
in Fig. 7. The panels show the median overdensity (top-left), tem- 
perature (top-right), hydrogen mass fraction (middle-left), metalli- 
city (middle-right), oxygen abundance (bottom-left), and O vi ion- 
isation fraction (bottom-right) as a function of fcovi and Nqvi- For 
reference, we have included contours showing the distribution (by 
number) of O vi absorbers in the fcovi - Nqvi plane. These contours 
contain, starting from the innermost, 25, 50, 75, 90, and 99 per cent 
of the identified O vi absorbers. 

Looking at the top-left panel, we can see that there is a 
mild correlation between column density and overdensity, with 
no apparent correlation between the line width and overdensity, 
which is consistent with our missing iovi - Nqvi correlation (see 
Fig. 2, bottom-right panel). In contrast, there is a strong correl- 
ation between Doppler parameter and the optical-depth weighted 
gas temperature, as shown by the top-right panel. Furthermore, the 
middle-right panel shows a good correlation between column dens- 
ity and gas metallicity, i.e. a metallicity-density relationship. This is 
consistent with the results shown in the middle-left and bottom-left 
panels, where we can see that the hydrogen mass fraction of the gas 
giving rise to O vi absorption decreases with Nqvi^ and correspond- 
ingly, the oxygen abundance increases. Finally, the bottom-right 
panel shows that there is no clear trend of the ionisation state of the 
O VI bearing gas with Z^ovi nor with Nqvi- 

All these results together imply that .some information about 
the physical conditions can be gained from the measured Ovi 
column densities and line widths alone, but only in a statistical 
sense, given the large scatter in the correlations between directly 
observable and derived physical quantities. In particular, the bot- 
tom panels of Fig. 7 reinforce our previous statement that assuming 
a fixed value for the oxygen abundance and O vi ionisation fraction 
to estimate the baryon content of O vi bearing gas is dangerous. 



5.2.1 Gas temperature, Doppler widths, and line strengths 

One of the main purposes of recent observational studies about in- 
tergalactic O vi absorption at low redshift has been to determine the 
temperature, and thus the ionisation state, of the gas giving rise to 
the observed absorption. Some of these (e.g. Thom & Chen 2008b; 
Tripp et al. 2008) have estimated the gas temperature using the line 
width of well-aligned, i.e., within some velocity uncertainty (t(Ai;), 
O VI - H I absorbers, and have found that the implied gas temperat- 
ure are typically T < 10^ K, indicating that the associated Ovi is 
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Figure 8. Comparison between iitted Doppler parameters and O vi gas tem- 
peratures for all components (3034) in sample 2. The filled dots show the 
median value in each bin of size Aio vi = 2 km s" ' (indicated by the A-bars), 
while the y-bars show the 25th and 75th percentiles, respectively. The blue 
dashed curve shows the relation between temperature and Doppler width 
assuming pure thermal broadening (see text for details). The (normalised) 
histogram in the sub-panel shows the distribution of foovi values. 



mainly photo-ionised. This approach, however, implicitly assumes 
that Ovi and Hi absorbers arise in the same gas, which might not 
be true in general, even for small velocity displacements. As also 
noted by Thom & Chen (2008b) and Tripp et al. (2008), the gas is 
most probable multi-phase, and a broader H i component related to 
the O VI absorption might be too weak to be detected, particularly 
if the metallicity of the component is high. Such a broad compon- 
ent would relax the upper limit on the temperature imposed by the 
Hi line width, leading to higher gas temperatures. Indeed, Tripp 
et al. (2008) have estimated an upper limit on the temperature of 
the O VI bearing gas using the O vi line widths alone and assuming 
pure thermal broadening, and they find significantly higher temper- 
atures than those allowed by well-aligned OVI - H i absorbers. 

Danforth & Shull (2008) find a strong correlation between the 
column densities of O vi and N v absorbers, and use this to model 
the ionisation state of the absorbing gas, concluding that the ob- 
served column density ratios are consistent with coUisionally ion- 
ised gas at r = 10^^*"' K, assuming a solar (N/0) abundance. 
Danforth & Shull (2008) note that although the observed Nqvi/Nnv 
ratios are also consistent with pure photo-ionisation models, the 
implied ionisation parameters, the required metallicities, and/or the 
spectral hardness of the photo-ionising radiation are not compatible 
with other constraints. 

Given the current disagreement about the temperature, and 
hence the ionisation state, of the gas traced by the observed O vi, we 
inspect the correlation between O vi line width and gas temperature 
suggested by the top-right panel of Fig. 7 in more detail. For this 
purpose, we bin the Tovi values using bins of size A^ovi = 2 km s"' 
for all identified O vi absorbers with foovi < 40 km s"' in sample 2, 
and compute the median, and 25th and 75th percentiles in each bin. 
The result is shown in Fig. 8. The relation between gas temperature 
and Doppler parameter assuming pure thermal broadening. 
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has been included in this plot for reference (blue dashed curve). 

Clearly, there is a tendency for the the optical-depth weighted 
gas temperature to increase with iovi- For feovi > 15kms"' (cor- 
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Figure 7. Physical conditions of tlie O vi bearing gas as a function of O vi observables, i.e. column density and Doppler parameter. The colour scale shows, 
for each (/:)ovi>Novi)-cell of size Afoovi = 2.0 kms"' and AlogNovi = 0.2 dex, the median overdensity (top-left), temperature (top-right), hydrogen mass 
fraction (middle-left), metallicity (middle-right), oxygen abundance (bottom-left), and O vi ionisation fraction (bottom-right). For reference, we have included 
contours showing the distribution (by number) of vi absorbers on the Z^ovi ~ Nqvi plane. These contours contain, starting from the innermost, 25, 50, 75, 
90, and 99 per cent of the identified O vi absorbers and they are the same in all panels. 



responding to log(r/K) > 5.3), non-thermal broadening becomes 
progressively more important, leading to larger i- values for a given 
temperature than expected for pure thermal broadening. Note that 
for the bin just above the resolution limit in our synthetic spectra, 
^min = 4.2 kms"', the median temperature is larger than expected 
for pure thermal broadening, as is the case for the 75th percentiles 
at higher fo-values, which is unphysical. As discussed in Sec. A2, 
this is an artefact of our fitting procedure in which weak, narrow 
components are added to broad, shallow absorption features to im- 
prove the fit. Some of these broader components arise in gas with 
relatively high temperatures and, correspondingly, high optical- 
depth weighted temperatures are assigned to the related (narrow) 
components. 

From the previous discussion we can conclude that line widths 
do reflect the temperature of the absorbing gas in our simulation, al- 
though in a statistical sense, but that no definite conclusion can be 
reached for individual absorbers, unless further information, e.g. 
other related absorbing ions, are available. Also, it is plausible that 
the correlation between gas temperature and line widths might be 



washed out by line-broadening mechanisms such as turbulence, 
which are not resolved by our simulation, or by a too low instru- 
mental resolution. 



Nevertheless, assuming we would estimate an upper limit on 
the gas temperature from the O vi line widths alone, e.g. by invert- 
ing the relation shown in Fig. 8, we would find a good agreement 
with the upper limits estimated by Tripp et al. (2008, their Table 7) 
from their O vi line widths, which are, however, much higher than 
the temperatures implied by their sample of well-aligned Ovi - 
Hi absorbers. This, in turn, suggests that gas temperatures derived 
from bnJhovi ratios alone might be too low. Furthermore, this sup- 
ports the idea that there are broad H i components which truly are 
associated with the O vi components, but which are too weak to be 
detected, as anticipated by Tripp et al. (2008). 



In addition, Ovi detections might be biased towards lower 
temperatures. This is plausible since absorption strength, as meas- 
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ured by the optical depth at the line centre'" , tq, is proportional to 
the ratio of the column density to the Doppler width, N jb. Hence, 
at a fixed column density, line strength is inversely proportional to 
the Doppler parameter and thus, assuming pure thermal broaden- 
ing, inversely proportional to the square-root of the gas temperat- 
ure. Also, at a fixed &-value, line strength is directly proportional to 
column density, and given the correlation between column density 
and metallicity presented in Sec. 5.1.2, also directly proportional 
to the gas metallicity. Since metallicity enhances cooling, it is thus 
natural to expect that stronger O vi absorption traces gas at lower 
temperatures. 

In order to test this, we bin the distribution of gas temperat- 
ures To VI as as function of the line depth (central optical depth) for 
our sample 2, and compute the median temperature and the 25th 
and 75th percentiles in each bin. The result is shown in Fig. 9. 
As we anticipated, we see a clear anti-correlation between the me- 
dian gas temperature and absorption strength. We note that Tripp 
et al. (2008) estimated their detection threshold to be Tq ~ 0.1. 
We can see from the histogram in the sub-panel - which shows the 
distribution of lines as a function of line strength - that there is 
a non-negligible number of weak O vi absorption features arising 
in warm-hot gas, which would not have been detected by Tripp 
et al. (2008). Qualitatively, this conclusion is neither affected by 
the assumed S/N (the same behaviour is found for spectra with 
S/N = 10), nor by the fact that in our synthetic spectra there is 
no apparent correlation between b and Nqvi- Incidentally, the lat- 
ter facilitates the interpretation of the anti-correlation between gas 
temperature and line strength shown in Fig. 9 as a signature of en- 
hanced cooling in enriched gas, which is discussed in the next sec- 
tion. The bottom line of this exercise is that observations, even if the 
gas temperature is measured accurately, may be biased towards ab- 
sorption systems tracing gas at somewhat lower temperatures than 
expected for collisionally ionised gas, due to the limited S/N of the 
data. 



5.3 Cooling times 

We presented in Sec. 5.1.1 our key result that simulated Ovi 
absorbers preferentially trace gas at temperatures in the range 
T ~ 10^ - 10^ K, and thus trace gas in the low temperature regime 
of the WHIM. This is in marked contrast with the results of Oppen- 
heimer & Dave (2009, their Fig. 14), who find that the majority of 
O VI absorbers arise in gas at temperatures lO'' * <T< lO'* ** K with 
a peak atT x lo"*-^*"- K. The purpose of this section is to address 
the source of this discrepancy. 

Oppenheimer & Dave (2009) find that in their simulations 
O VI traces over-enriched (by factors of four to six) regions with 
a clumpy metallicity distribution, and they argue that these re- 
gions are thus subject to enhanced metal-line cooling, such that 
Ovi bearing gas which was initially shock-heated to temperat- 
ures T > 10* K is able to cool to photo-ionised temperatures 
well within a Hubble time. Our results are consistent with Oppen- 
heimer & Dave (2009) inasmuch as we also find that Ovi traces 



'" The central optical depth of an absorption line with column density N 
and Doppler parameter b is given by 



Nqvi *ovi [cm km s] 
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where e and nig are the electron's charge and mass, A and / are the trans- 
ition's rest wavelength and oscillator strength, respectively, and the last nu- 
merical factor is valid for the Ovi /tl031.93 A transition. 




Line Depth (tq) 

Figure 9. Optical depth-weighted gas temperature of the O vi absorbing 
gas as a function of central optical depth tq (or line strength Nqvi/^Ovi) 
for all components identitied in 5000 synthetic spectra with S/N = 50 at 
z = 0.25. The blue dashed lines indicate the median gas temperature at the 
smallest and largest line depths. The (normalised) histogram in the lower 
panel shows the distribution of O vi components contributing to each bin. 
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Figure 10. Distribution of O vi mass (coloured regions) and distribution of 
O VI absorbers (white contours) in the temperature-overdensity plane. Note 
that these distributions are identical to those shown in the bottom panel of 
Fig. 3. Cooling time contours (magenta lines) for two different metallicities 
(solar, solid; 10 per cent solar, dashed) split the plane into regions con'es- 
ponding to temperatures and overdensities for which the cooling time is 
shorter (upper part) and longer (lower part) than a Hubble time. 



in-homogeneously enriched gas with relatively high metallicities, 
but at significantly higher temperatures. In other words, our results 
indicate that, in our simulation, the vast majority of the gas traced 
by Ovi that was shock-heated to temperatures T > 10^ K has not 
yet cooled down to photo-ionisation temperatures, in spite of its 
enhanced metallicity. 

This fact is demonstrated in Fig. 10, where we reproduce the 
bottom panel of Fig. 3 showing the O vi mass distribution (colour 
shading) and the distribution of O vi absorbers (white contours) in 
the r - A phase diagram, together with the locus of (T, A)-values 
for which the cooling time equals the Hubble time for gas of solar 
(solid) and ten per cent solar (dashed) metallicity. Note that the 
metallicity values used roughly span the range of gas metallicities 
traced by our Ovi absorbers (see Fig. 4). Here, the cooling time is 
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Figure 11. Normalised, absolute net cooling rates at z = 0.25 for gas with nn = 10" cm" , for gas with metallicity of 10 per cent solar (left panel) and solar 
(right panel), computed assuming collisional ionisation equilibrium (CIE; red), including photoionisation (Wiersma et al. 2009a, PI; black), and including 
photoionisation only for gas of primordial composition while assuming CIE for metal-line cooling (Hybrid; blue dashed). See text for details. 



defined as 



'cool 



dT/dt 



3 nkT 



(5) 



where A„et is the net {i.e., cooling minus heating) radiative cool- 
ing rate at a fixed density. Note that the contour segments to the 
left from T » 10"*^ K correspond actually to heating times. The 
point where heating and cooling time contours merge defines the 
equilibrium temperature, i.e. the temperature at which the net cool- 
ing rate vanishes and the cooling time becomes effectively infinite 
(see also Fig. 11). From Fig. 10, it becomes apparent that, while 
virtually all O vi absorbers in our simulation at z = 0.25 will cool 
down to the equilibrium temperature in less than a Hubble time" 
(in the absence of further shock-heating), most will still be hot by 
z = Q. Note that, as anticipated in Sec. 5.1.1, the low-temperature 
O VI bearing gas has nearly reached the equilibrium temperature 
T ~ 10"" K, typical for photo-ionised O vi. 

The above result strongly suggests that an important factor 
leading to the notable difference between the results of Oppen- 
heimer & Dave (2009) and our findings regarding the mean tem- 
perature of gas traced by O vi is the way the radiative cooling of gas 
is treated in the simulations. As described in Oppenheimer & Dave 
(2009), metal-line cooling is included in their simulations using the 
models by Sutherland & Dopita (1993), which assume collisional 
ionisation equilibrium and fixed relative abundances. They only in- 
cluded photo-ionisation for hydrogen and helium. In contrast, as 
described in detail in Schaye et al. (2010), the OWLS runs include 
radiative cooling implemented according to the method by Wi- 
ersma et al. (2009a), who compute cooling rates on an element-by- 
element basis in the presence of the cosmic microwave background 
(CMB) and an ionising (i.e., UV/X-Ray) radiation field as modelled 
by Haardt &. Madau (1996). Wiersma et al. (2009a) showed that 
including photo-ionisation by the meta-galactic UV/X-Ray back- 
ground not only drastically decreases the cooling rates for hydrogen 
and helium, but also for heavy elements. Moreover, they showed 
that photo-ionisation of heavy elements shifts the equilibrium tem- 
peratures to higher values. In other words, photo-ionisation can 



For reference, the Hubble time thuI 



^(z = 0) = 13.8 Gyr and 



c(z = 0.25) = 10.9 Gyr for our adopted cosmology. 



drastically affect the radiative cooling rates of enriched gas and cal- 
culations neglecting photo-ionisation will thus underestimate the 
temperature of the gas. 

A comparison of the cooling rates computed according to each 
method is shown in Fig. 11. In each case, the gas is exposed to 
the CMB and the UV/X-Ray background at z = 0.25, and has a 
hydrogen density Uh = 10"^ cm"' (A « 50), a hydrogen mass frac- 
tion Xh = 0.752 - which roughly correspond to the median values 
of (nH)ovi and (Xh)ovi in sample 2 - and a metallicity of 10 per 
cent solar (left panel) and solar (right panel). The black curve in 
either panel shows the net cooling rate including photo-ionisation 
of heavy elements (PI); the red curve, the net cooling rate assum- 
ing pure collisional ionisation equilibrium (CIE) and fixed (solar) 
abundances, while the blue curve shows the net cooling rate as- 
suming CIE, but including photo-ionisation for hydrogen and he- 
lium (Hybrid). The dip apparent in each cooling curve shows the 
equilibrium temperature. Note, in particular, the difference in equi- 
librium temperatures and overall amplitude between the black and 
blue curves, which correspond to the cooling rates implemented in 
the OWLS runs and the approach adopted by e.g., Oppenheimer & 
Dave (2009), respectively. It is clear that simulations using the lat- 
ter prescription overestimate the cooling rates. The difference in the 
overall amplitude and in the equilibrium temperature is more severe 
at lower densities ~ over an order of magnitude at nn = 10"*cm"' 
- and higher metallicities (see right panel). 

The reduction in the cooling rates and the corresponding in- 
crease in the cooling times, together with the shift of the equi- 
librium temperature to higher values, may explain the difference 
between the results of Oppenheimer & Dave (2009) and ours re- 
garding the temperature of gas traced by simulated O vi absorbers. 
This difference is even more striking considering the fact that Op- 
penheimer & Dave (2009) use particle abundances while we use 
smoothed metallicities throughout the simulation to compute the 
metal cooling rates. As previously mentioned and discussed in de- 
tail by Wiersma et al. (2009b), the use of smoothed abundance res- 
ults in increased cooling for most of the enriched gas, although the 
cooling rates are reduced for gas particles that have metal mass 
fractions higher than those of their neighbours. One would thus 
naively expect to find exactly the opposite result, namely that our 
simulated O vi absorbers trace gas at lower temperatures. That this 
is not the case indicates that the reduction of the radiative cooling 
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rates due to photo-ionisation has a large impact on the thermal state 
of the gas. 

One should keep in mind that our fiducial run is not directly 
comparable to the set of simulations by Oppenheimer & Dave 
(2009), since they differ in various aspects, in particular in the 
assumed cosmology (WMAP-5 vs. our assumed WMAP-3), and 
in the sub-grid prescriptions to model, e.g., SN feedback in the 
form of galactic winds. Also, recall that we use different meth- 
ods to ascribe physical properties to absorption features identified 
in synthetic spectra. However, we have carried out runs in which 
we have varied the cosmology, the prescriptions for SN feedback, 
and even have implemented the method by Oppenheimer & Dave 
(2009) to estimate the gas temperature and density related to O vi 
absorption features, and our results remain practically unaltered. 
We will report on such variations to our reference model in a future 
paper where we analyse the physical conditions of O vi absorbers 
in several of the OWLS runs, and we will show that the fact that 
Ovi traces gas at temperatures r~10^-10*K and overdensit- 
ies 1 « A < 10' is a robust prediction of these simulations. This 
demonstrates the capital importance of a correct treatment of radi- 
ative cooling in simulations that attempt to shed light on the much 
debated physical state of absorbing gas giving rise to absorption by 
highly ionised species, currently observed in QSO/AGN spectra. 
We should therefore keep in mind that, although our implementa- 
tion of radiative cooling is more advanced than the ones used before 
in that we include photo-ionisation of heavy elements by the UV 
background and in that we do not assume fixed relative abundances 
of different elements, we are still assuming ionisation equilibrium 
and we are still neglecting local sources of ionising radiation. 



6 CONCLUSIONS AND DISCUSSION 

It is now well established that the inventory of observed baryons 
in the low-redshift Universe falls short of the baryonic mass meas- 
ured at high redshift and expected from big bang nucleosynthesis 
calculations. However, numerical simulations of structure forma- 
tion predict that at low redshift a significant amount of baryons 
are hidden in a hard-to-observe, diffuse, hot gas phase termed the 
Warm-Hot Intergalactic Medium (WHIM), thus providing a poten- 
tial solution for the mismatch between observed and expected ba- 
ryons. As a consequence, a large observational campaign has been 
launched with the primary goal of detecting this gas phase and es- 
timating its baryon content, mainly through the observation of O vi 
seen in absorption in the spectra of QSO and AGN. 

Although statistically significant samples of low-redshift O vi 
absorption systems have been obtained by way of UV absorption 
spectroscopy, the nature and the physical state of these absorbers 
remains controversial. Due to the complexity of the absorbing gas 
in terms of its density, temperature, and ionisation state, no secure 
estimate of the amount of baryons contained in this gas has been 
achieved. Nor is it clear in which gaseous phase these absorbers 
arise and whether they are predominantly collisionally or photo- 
ionised. Observations thus need to be supported by theoretical work 
and, given the intricate interplay between the various physical pro- 
cesses governing the formation and evolution of metal absorption 
systems, the use of numerical simulations has become a powerful 
tool to complement (and/or reinforce) the results inferred from ob- 
servational data. 

In this paper, we have made use of a set of cosmological hy- 
drodynamical simulations of the 'reference model' from the Over- 
whelmingly Large Simulations (OWLS) project (Schaye et al. 



2010) to study the physical properties and baryon content of in- 
tervening Ovi absorbers at low redshift. The main advantage of 
these simulations compared to previous ones is the implementation 
of important physical processes that have been largely ignored in 
earlier studies. The delayed release of mass by massive stars, inter- 
mediate mass stars, and supernovae of types la and II is followed 
for all the 1 1 elements that contribute significantly to the radiative 
cooling rates. Radiative cooling is computed element-by-element 
and the effect of photo-ionisation by the UV/X-ray background is 
taken into account not only for hydrogen and helium, but also for 
the heavy elements. Wiersma et al. (2009a) showed that the sup- 
pression of cooling by photo-ionisation of heavy elements, which 
was ignored in previous studies, is important for the WHIM and we 
confirmed that it has large consequences for the gas that is respons- 
ible for the O vi absorption. 

We have generated a large number of synthetic spectra that 
closely match the properties (S/N, resolution, etc.) of observed 
spectra. We have identified and fitted individual Ovi absorption 
components, and have measured their column density, Doppler 
parameter, and rest equivalent width, and their corresponding dis- 
tributions. We have shown that we can reliably identify lines with 
O VI column densities above a given threshold that increases with 
decreasing signal-to-noise ratio, even though we still miss a frac- 
tion of the (weak) absorption lines present in the spectra. Moreover, 
we showed that we are able to measure rest equivalent widths quite 
accurately, and that the fitted Doppler parameters are reliable for 
moderate signal-to-noise ratios. Finally, we showed that the line 
parameter distributions are converged with respect to box size and 
resolution. 

Our simulations reproduce the shape of the observed column 
density distribution function, although with a lower (by a factor 
~ 2) amplitude. Similarly, the predicted cumulative equivalent 
width distribution has a lower amplitude, a lack of strong lines, and 
an excess of low equivalent width lines. Furthermore, our simula- 
tions do not show the positive correlation between fcovi and Nqvi 
that appears to be observed for O vi absorbing systems in a variety 
of environments, although the observations are uncertain and even 
contradictory in this respect. We attribute the lack of strong (i.e., 
high equivalent width) lines in our simulation to a lack of lines 
with high O vi column density and large Doppler parameter values. 

A similar disagreement between the predicted and the ob- 
served Ovi cumulative line-number density and fcovi ~ Nqvi cor- 
relation has been reported by Oppenheimer & Dave (2009), who 
argued that small-scale turbulence is the crucial mechanism leading 
to the large fc- values measured in IGM O vi systems, and solved the 
discrepancy between their simulations and observations by adding 
sub-resolution turbulent broadening to the Doppler parameters of 
their simulated absorption lines. Exploring the effects of such a 
post-run variation on the physical properties of the absorbing gas 
is beyond the scope of the present work, but we will do so in a 
future paper 

We stated that even with simulations it is not trivial to as- 
sociate physical properties, such as a density, temperature, chem- 
ical composition and ionisation balance, to a particular absorber 
because gas parcels with a range of properties will generally con- 
tribute to a single absorption line. Therefore, we chose to assign 
physical properties to each pixel in velocity space by weighting all 
gas parcels by their contributions to the O vi optical depth in that 
pixel. The resulting optical-depth weighted quantities were then in 
turn averaged over the line profile, weighted by the optical depth 
of each pixel. The use of optical-depth weighted quantities allows 
the direct association of an absorption feature with the gas where 
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the absorption takes place, even in the presence of redshift-space 
distortions resulting from peculiar velocity gradients and thermal 
broadening. 

In terms of the physical state and gas (metal) mass content of 
the O VI bearing gas, we have found that: 

(i) Detectable Ovi absorbers trace gas at temperatures 
W-^ K < r < 10" K with a peak at T « 10'-' K, and overdensities 
1 « A < 10^ 

(ii) The bulk of the O vi mass in our simulation is found in gas at 
moderate overdensities A ~ 10"^ - 10^^, and temperatures between 
r = 10"* K and 10* K. While approximately 30 per cent of the O vi 
seen in absorption arises in gas at T < 10^ K, the rest is found to 
trace a the hotter gaseous phase at temperatures around T x 10'-' K. 

(iii) The metallicity distribution of the gas traced by Ovi is 
narrow, extending over the range (0.1, 1) Zg, with a peak at 
Zovi ~ 0.6 Zq, i.e., almost an order of magnitude higher than the 
median metallicity of gas at comparable overdensities and temper- 
atures and also about an order of magnitude higher than would be 
inferred from the column density ratio of O vi to H i for aligned 
O vi-H I pairs. 

(iv) The Ovi bearing gas contains a baryon fraction 
Q.h(Ovi)/iih x: 0.014, and contains less than one per cent of 
the total metal mass in our simulation. 

The low baryon fraction associated with the Ovi absorbers 
is mainly a consequence of the fact that the metal distribution is 
patchy. The Ovi absorption arises mostly in high-metallicity re- 
gions that are small compared with the scales on which the dens- 
ity and the temperature fluctuate. Hence, if we were to recompute 
the synthetic spectra after smoothing the metal distribution, we 
would probably find that O vi traces a larger fraction of the bary- 
ons. The small-scale inhomogeneity of the metal distribution there- 
fore makes the question of which baryons to associate with the O vi 
gas somewhat ill defined. This is particularly relevant because the 
lack of small-scale metal mixing may be partly numerical. How- 
ever, there is in fact strong observational evidence that the metals 
in the real Universe are also poorly mixed on small scales (Schaye 
et al. 2007) and that the metal distribution is indeed patchy on the 
scales over which the H i absorption is uniform. 

The common approach of comparing Ovi and Hi column 
densities to estimate the physical conditions in intervening ab- 
sorbers from QSO observations may thus be misleading, as most 
of the H I (and most of the gas) is not physically connected with 
the high-metallicity patches that give rise to the O vi absorption. A 
possible solution to this problem is to rely only on upper limits for 
the absorption by H i (and other species) at the velocity of O vi, as 
was done for Civ absorbers in Schaye et al. (2007). 

The physical conditions in the absorbing gas can in prin- 
ciple be inferred from the measured O vi column densities and line 
widths, but only in a statistical sense, given the large scatter in 
the correlations between directly observable and derived physical 
quantities. In particular, the gas temperature can be recovered reas- 
onably well from the O vi line widths for fo-parameters < 15 km s"' 
(log(r/K) < 5.3), but for larger line widths non-thermal broaden- 
ing becomes increasingly important. 

The conclusion that Ovi absorption preferentially traces gas 
at r ~ 105.3+0.5 j^ typical of coUisionally ionised gas, is in marked 
contrast to the results of Oppenheimer & Dave (2009), who find 
that the overwhelming majority of the O vi absorbers in their sim- 
ulations trace gas at temperatures T « j (34.2+0.2 j^^ which are typ- 
ical of photo-ionised gas. We showed that this drastic difference 
is most probably a consequence of the different treatment of ra- 



diative cooling of enriched gas. In particular, while Oppenheimer 
& Dave (2009) assumed collisional ionisation equilibrium for the 
heavy elements when computing cooling rates, we included the ef- 
fect of ionisation by the UV/X-ray background radiation, which 
strongly suppresses the cooling rates for low-density gas. 

While the mismatch between our predicted and observed line 
parameter distributions is certainly unsatisfying, we think that it 
may not undermine our conclusions about the physical properties 
of the O VI absorbing gas. As Oppenheimer & Dave (2009) showed, 
it is possible to alleviate the discrepancy with observations by using 
a physically motivated post-run modification, namely the addition 
of turbulence on scales that are unresolved by the simulation. We 
do not expect such a modification to significantly affect the phys- 
ical properties derived for the absorbing gas, precisely because it is 
done in a post-run fashion. 

It should also be stressed that it would a priori be surprising 
if we matched the column density distributions very well, given 
that the nucleosynthetic yields are uncertain at the factor of a 
two level, even for a fixed stellar initial mass function (e.g. Wi- 
ersma et al. 2009b). Some rescaling of the abundances in post- 
processing would thus be justified, although this would break the 
self-consistency of the simulation if it were to change the cooling 
rates significantly. 

We postpone a detailed presentation and discussion of a series 
of modifications to our reference run (i.e., simulation runs with dif- 
ferent parameters) as well as post-run variations to a future paper. 

The results presented in this work clearly show that O vi ab- 
sorbers trace neither the main metal reservoirs nor the bulk of the 
baryons in our simulation. Instead, they trace over-enriched, shock- 
heated gas at moderate overdensities. Hence, our study indicates 
that Ovi - albeit coUisionally ionised - is a poor tracer for the 
bulk of the baryons in the WHIM, as most of the baryons reside 
in shock-heated gas that has much lower metallicities. Instead, our 
simulations indicate that O vi absorbers may trace cooling metal- 
ejecta from galaxies and can thus be used to study the circulation 
of metal-enriched gas in the intergalactic environment of galaxies. 
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APPENDIX A: ABSORPTION LINE PARAMETERS 

In this section we explore tiie effect of the signal-to-noise ratio 
(S/N) on the distributions of O vi line parameters (column density, 
Doppler parameter, rest equivalent width) obtained from the ana- 
lysis of 5000 synthetic spectra in the range z = 0.0 - 0.5 which 
span a total redshift path Az = 98.5 (sample 1; see Sec. 3.1). We 
analyse each of these distributions in detail to assess the strengths 
and weaknesses of our spectral analysis. 



A I Column density distribution 

We first consider the distribution of Ovi column densities, ob- 
tained by binning the set of fitted O vi column densities for each 
adopted S/N in logarithmic bins of size AlogNovi = 0.2 dex, 
and normalising by the total redshift path surveyed. In each bin, 
we compute Poisson single-sided Icr confidence limits using the 
tables by Gehrels (1986). Note that this uncertainty does not in- 
clude the fitting error in the column density, which is typically 
around 10 per cent. The histograms thus obtained are displayed 
in the top panel of Fig. Al. This plot clearly shows the effect of 
the S/N on the detection limit in terms of column density: First, a 
decreasing S/N leads to a smaller number of identified components 
(3115 at S/N = 50 vs. 1093 at S/N = 10). Second, it increases the 
column density threshold - log(Novi/cm"^) » 12.3 at S/N = 50 
vs. log(Novi/cm"^) K 13.2 at S/N = 10 - above which lines can 
be reliably detected. The column density detection threshold for 
each S/N can be best appreciated in the inset displaying the norm- 
alised cumulative distribution, i.e., the number of lines with column 
densities above a given value. Note that a column density threshold 
log(Novi/cm"^) a 13.2 for S/N = 10 is typical for Ovi absorbers 
identified in FUSE and HST/STIS data at a similar S/N (e.g. Tripp 
et al. 2000; Danforth & ShuU 2005; Thom & Chen 2008b). 

In order to assess the completeness of our line sample and 
to check the quality of our fits, we compare the integrated Ovi 
column density along each physical sightline in the simulation (the 
true total O vi column density) to the sum of the column densit- 
ies of all O VI absorption components identified in the correspond- 
ing spectrum (the recovered total O vi column density). The res- 
ult is shown in the bottom panel of Fig. Al. We bin the relat- 
ive difference (5Novi/Novi of each sightline (spectrum) as a func- 
tion of the true column density using A log Nqvi = 0.2 dex (as dis- 
played by the x error-bars), and compute the 25th and 75th percent- 
ile in each bin (corresponding to the lower and upper y error-bar, 
respectively). Comparison to a perfect match given by the black 
dashed line shows that the agreement between the true and the re- 
covered column densities is good, although with a significant scat- 
ter, and with the true amount of Ovi along each physical sight- 
line being underestimated at all column densities but the lowest. 
In particular, we can see that we miss ~ 30 per cent of the O vi 
column density along sightlines with a total column density around 
No VI = 10'^^ cm"-. This is due to the fact that along some of these 
sightlines (~ 20 per cent), more than one absorption component 
(usually two or three) are identified. Hence, some of the absorption 
components along these sightlines have column densities that fall 
below the detection threshold. Of course, this is true in general for 
all sightlines, which explains why we systematically underestim- 
ate the total O vi column density. Note that the particular column 
density at which the recovery is worst depends on the adopted S/N. 
Moreover, at high column densities, saturation effects may also 
play a role. Indeed, since Ovi saturates at Nqvi ~ 2x 10 '"'cm"- 
for iiQvi = 4.2 kms"' {i.e., at our chosen spectral resolution limit; 
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Figure Al. Top Differential column-density distribution 
c9^A/'/(51ogNovi&) (histograms) for different S/N values. The error 
bars show Poisson single-sided Icr confidence limits, computed according 
to the tables by Gehi'els (1986). The inset displays the corresponding 
normalised cumulative distributions, i.e., the number of lines above a 
given column density. Bottom: Comparison between the integrated column 
density along each physical sightline (the true column density) and the 
sum of the column densities of the fitted lines in the corresponding 
spectrum (the recovered column density). The y-axis displays the relative 
dift'erence between true and recovered column densities, (5Novi / Nqvi- 
We have binned the results using A log No vi = 0.2 dex (as indicated by the 
X error-bars); the dots and y error-bars show the median value in each bin, 
and 25 and 75th percentiles, respectively. The black dashed line shows a 
one-to-one con'espondence. Colours as in the top panel (see legend). 



see next section), it is possible that lines with true column densities 
above this value are fitted using components with somewhat lower 
column densities. The latter effect is, however, negligible, given the 
small number of components with Nqvi > 2 x 10 '* cm"^ (see inset 
in top panel of Fig.Al). 

A2 Doppler parameter distribution 

Since the STIS resolution limit of FWHM = 7kms"' translates 
into a minimum Doppler parameter /7,^jn = 4.2 km s"', we impose 
the condition that our fitting algorithm discards any identified fea- 
tures with fitted i- values smaller than this to constrain the detection 
of spurious features. However, this does not have a large impact 
on our resulting line width statistics, given the small value of the 
threshold relative to the typical O vi line widths (Z7ovi ~ 12 km s"' ; 
see below). Note that for other instruments with lower resolution 
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Figure A2. Distribution of Doppler parameters (histograms). Note tliat the 
distributions have been normalized to unit area. The inset shows the norm- 
alised cumulative distribution for each S/N. The black, dashed line shows 
the STIS resolution limit at FWHM = 7 kin s"' , which translates into a min- 
imum Doppler parameter &,„!„ = 4.2 kin s" ' . The cut-off at this value is due 
to the condition imposed on our fitting algorithm that any line with fitted 
fc-value smaller than fomin be discarded. 



such as COS (FWHM « 17 km s"'), imposing a minimum allowed 
fe-value in the fitting procedure might have a noticeable effect on 
the line width distribution. 

The distributions of Doppler parameters obtained from our 
spectra for S/N = 10, 30, and 50 are shown in Fig. A2. As can be 
appreciated, all three distributions are similar in width, with peaks 
at Doppler values in the range foovi ~ 10 - 12kms"', and all show 
an extended tail towards fc-values as large as 50 kms"'. The me- 
dian Doppler parameter is i,ned = 12.9, 12.7, and 11.7 kms"' for 
S/N = 50, 30, and 10, respectively, and typical relative errors in the 
fitted Doppler parameters are around 10 per cent, independent of 
the adopted S/N. The error bars display Poisson uncertainties. We 
find that a low S/N leads to a slight decrease in the number of com- 
ponents with feovi > 10 km s"', as can be seen in the inset display- 
ing the normalised cumulative distribution for each S/N. This is 
expected, since broader components are generally shallower (for 
a fixed equivalent width) and get progressively lost in the noise 
with decreasing S/N. To a lesser extent, it is also true that with 
decreasing S/N, intrinsically broad lines are occasionally contam- 
inated by strong noise spikes and may be thus be fitted with mul- 
tiple narrow components, thus effectively improving the fit. This 
can have important (and unwanted) consequences for the correla- 
tion between the temperature assigned to an absorption line [i.e., 
to the gas where the absorption takes place) and its corresponding 
Doppler parameter. Note on the other hand that, irrespective of the 
S/N, a high fraction (more than 70 per cent) of the fitted compon- 
ents have Doppler widths 10 km s"' < fcovi < 40 km s"' , which for 
oxygen coiTespond to temperatures 10'' K < T < 10* K (assuming 
pure thermal broadening; see equation 4). The relation between 
Doppler parameter and gas temperature is discussed in detail in 
Sec. 5.2.1. The similarity between the Doppler parameter distribu- 
tions for S/N = 30 and S/N = 50 suggests that our fitted Doppler 
parameter are reliable for S/N > 30. 



A3 Equivalent width distribution 

Here we further explore the quality of our fits by analysing the 
measured rest equivalent widths W^. The rest equivalent width of 



an O VI absorption line is measured directly by integrating the flux 
depression over the entire spectrum, and correcting for the line 
broadening due to cosmological expansion. In the top-left panel 
of Fig. A3, we compare the rest equivalent width as a function of 
column density for each O vi component identified in spectra with 
different S/N to theoretical curves-of-growth (CoG) for a range of 
Doppler parameters which roughly corresponds to the width of the 
distributions shown in Fig. A2. It is remarkable that the overwhelm- 
ingly majority of the identified components lie in the range allowed 
by the curves-of-growth. The exception are a few absorption lines 
for which either the measured equivalent width is systematically 
underestimated (assuming the column density has been accurately 
measured), or the column density is slightly and systematically 
overestimated. The good correspondence between rest equivalent 
width and column density as compared to the CoG is even more 
notable if one considers that the rest equivalent width of a given 
line has not been determined from a CoG analysis, i.e., its meas- 
urement is completely independent of the corresponding fitted O vi 
column density and the fitted &-value. This demonstrates further the 
accuracy of the line fitting algorithm we use. We note that most of 
the absorption components, regardless of the S/N, lie on the linear 
part of the CoG. For these lines, the rest equivalent width scales 
linearly with the column density, and is (nearly) independent of the 
Doppler parameter. Unfortunately, this degeneracy makes it diffi- 
cult to assess the goodness of the Doppler parameter fitted to a par- 
ticular line. On the other hand, this explains why all identified O vi 
components seem to have (Nqvi. W,)-values consistent with CoG 
for foovi = 4-50 km s"' , even though we do find a few lines with 
(possibly artificial) ft- values as high as iiovi = 60 km s"' . 

Given the measured rest equivalent widths W ,, we compute 
the cumulative equivalent width distribution, i.e., the frequency of 
absorption lines per unit redshift above a given equivalent width 
threshold, for S/N = 10, 30, and 50. As already discussed in Sec. 4, 
we explicitly distinguish between individual absorption compon- 
ents and systems in order to have a better comparison to different 
sets of observations. The cumulative equivalent width distributions 
for individual components and systems adopting different S/N are 
shown in the top-right panel of Fig. A3. A lower S/N leads to a 
smaller number of lines at rest equivalent widths Wj < 30mA. 
The reason for this is two-fold: First, weak absorption features 
simply get lost in the noise. Second, with decreasing S/N, otherwise 
resolved narrow components are occasionally fitted by a single, 
broader absorption feature. Both these effects reduce the relative 
number of components at low equivalent widths and increase their 
number at high equivalent widths. The cumulative equivalent width 
distribution for systems shows qualitative trends which are similar 
to those of the distribution of individual components at low equi- 
valent widths. However, an important difference is that the number 
of identified systems seems to be nearly independent of the adop- 
ted S/N at the high equivalent-width end. This can be understood in 
terms of the definition of systems together with the fact that noise 
has little effect on strong absorption features. Note also that there 
are more strong {i.e., at Wr ^ 50mA) systems than strong compon- 
ents, as expected, since the latter are a subset of the former. 

While a cumulative distribution gives an idea about the amp- 
litude of the overall equivalent width distribution, a differential dis- 
tribution gives insight into the shape of the distribution. It is thus of 
interest to take a look at the differential distribution of rest equival- 
ent widths both of components and systems at different S/N values. 
The bottom-left panel of Fig. A3 shows the distribution of compon- 
ents (normalised to unit area). Here we see again that low W , lines 
get systematically lost with decreasing S/N, while the opposite is 
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Figure A3. Top left: Rest equivalent width Wi as a function of O vi column density for each O vi component and for different S/N (coloured symbols). The solid 
curves show theoretical curves-of-growth (CoG); the fe-value used to compute each CoG is indicated along the corresponding curve. Top right: Cumulative 
equivalent width distribution, i.e., number of lines per unit redshift above a given rest equivalent width for components (solid) and .systems (dashed) identified 
in our synthetic spectra for dift'erent S/N. Bottom left: Frequency of absoiption lines per equivalent width bin for individual components and dift'erent S/N. 
EiTor bars denote Poisson uncertainties. Bottom right: Same as bottom-left panel, but for absorption systems. 



true for high W, hnes, owing to the fact that individual features 
are occasionally merged in noisy spectra. This suggests that some 
of the high W, hnes identified in real QSO spectra could be artifi- 
cial. The bottom-right panel shows the differential W, distribution 
of systems. Again, weak systems seem to get progressively lost in 
the noise'^, as is the case for single components, even though the 
efi^ect is almost not present for S/N > 30, as can be judged by com- 
paring the distributions for S/N = 30 and S/N = 50 which are nearly 
indistinguishable. This indicates that, above a certain S/N value, the 
identification of absorption systems is less sensitive to the adopted 
S/N as compared to single components. 



APPENDIX B: CONVERGENCE TESTS 

In this section we address the numerical convergence of our results 
with respect to box size and resolution, focusing on the predicted 
O VI line parameter (column density, Doppler parameter, differen- 
tial equivalent width) distributions. To this end, we generate and 
analyse 5000 synthetic spectra at z = 0.25, adopting a standard S/N 



'- Recall that the detection threshold To, used to identify absorption sys- 
tems depends on S/N. 



of 50, for different runs including our reference model. Note that 
similar convergence tests in terms of the predicted star-formation 
history (SFH) and cosmic metal distribution have previously been 
presented and discussed in detail by Schaye et al. (2010, their sec- 
tion 3.2) and Wiersma et al. (2009b, their appendices B and C), 
respectively, using the same simulation runs that we use here. For 
this reason, we will restrict ourselves to the presentation of our con- 
clusions and kindly refer the reader to the studies mentioned above 
for a thorough discussion. 

Bl Box size 

Fig. B 1 shows the O vi column density distribution function (top- 
left), Doppler parameter distribution (top-right), and differential 
equivalent width distribution for single components (bottom-left) 
and systems (bottom-right), obtained from spectra synthesised us- 
ing three different simulations: REFJD25N128, REFJ.050N256, 
and REF_L100N512. All these simulations have the same spatial 
and mass resolution (see Table 1), but the sizes of the simulation 
boxes differ by factors of two. This plot shows that the prediction 
from our reference run for the O vi line parameter distributions is 
robust with respect to the size of the simulation volume. A box of 
L = 25 h^' Mpc on a side is sufficiently large to compute a reliable 
prediction for all three observables for lines with W^ < 200 mA. 
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Figure Bl. Convergence of our reference model with respect to simulation box size (at fixed spatial and mass resolution; see Table 1 for details). Top- 
left: O VI column density distribution. Top-right: O vi Doppler parameter distribution. Bottom-left: O vi equivalent width distribution for single components. 
Bottom-right: O vi equivalent width distribution for systems. 

Note, however, that we requke a box of at least 50 hr^ Mpc to cor- 
rectly sample the rarer strong {i.e., W,. > 500 mA) absorbers. 

B2 Resolution 

We now explore the effect of varying resolution on the pre- 
dicted Ovi observables. For this purpose. Fig. B2 compares a 
low resolution run REFJJ)50N128 and a moderate resolution 
run REF_L050N256 to our fiducial high-resolution simulation 
REFJJ)50N512. Note that we keep the box size fixed in this com- 
parison. While the low-resolution run is clearly different, the in- 
termediate resolution model is nearly identical to our fiducial run, 
suggesting that the results are converged with respect to the nu- 
merical resolution, except perhaps for the strongest absorbers with 
W, > 200 mA. 
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Figure B2. Convergence of our reference model (and fiducial run) with respect to resolution (at fixed box size; see Table 1 for details). Top-left: O vi column 
density distribution function. Top-right: O vi Doppler parameter distribution. Bottom-left: O vi equivalent width distribution for single components. Bottom- 
right: O VI equivalent width distribution for systems. 
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